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Topological order in two dimensions can be described in terms of deconfined quasiparticle exci¬ 
tations - anyons - and their braiding statistics. However, it has recently been realized that this 
data does not completely describe the situation in the presence of an unbroken global symmetry. 
In this case, there can be multiple distinct quantum phases with the same anyons and statistics, 
but with different patterns of symmetry fractionalization - termed symmetry enriched topological 
(SET) order. When the global symmetry group G , which we take to be discrete, does not change 
topological superselection sectors - i.e. does not change one type of anyon into a different type 
of anyon - one can imagine a local version of the action of G around each anyon. This leads to 
projective representations and a group cohomology description of symmetry fractionalization, with 
H 2 (G,A) being the relevant group. In this paper, we treat the general case of a symmetry group 
G possibly permuting anyon types. We show that despite the lack of a local action of G , one can 
still make sense of a so-called twisted group cohomology description of symmetry fractionalization, 
and show how this data is encoded in the associativity of fusion rules of the extrinsic ‘twist’ de¬ 
fects of the symmetry. Furthermore, building on work of HermelJ^, we construct a wide class of 
exactly solved models which exhibit this twisted symmetry fractionalization, and connect them to 
our formal framework. 
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I. INTRODUCTION 

In the last 30 years it has been realized that there 
exist quantum phases of matter that - unlike ordinary 
crystals, magnets, or superconductors - cannot be un¬ 
derstood in terms of symmetry breaking and local order 
parameters. The main example is the fractional quan¬ 
tum Hall effect, which instead exhibits a subtle non-local 
order manifested in exotic properties like emergent ex¬ 
citations with exotic statistics (anyons), protected gap¬ 
less edge modes, and ground state degeneracy on sur¬ 
faces of non-trivial topology. In general, one can use the 
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braiding statistics of the anyons to give a coarse classi¬ 
fication of gapped Hamiltonians; this is called intrinsic 
topological order. Although intrinsic topological order 
is independent of any symmetry considerations, it has 
recently been realized that the presence of a global sym¬ 
metry can further refine the coarse classification given by 
intrinsic topological order. In particular, there can exist 
several ‘symmetry protected’ quantum phases (SPTs) re¬ 
alizing trivial intrinsic topological order 2 "' 5 , and multiple 
‘symmetry enriched’ phases (SETs) corresponding to the 
same intrinsic topological order, with the latter being 
the focus of this paper. Our main result concerns an un¬ 
conventional type of symmetry action where acting with 
the global symmetry on some excitations turns them into 
new excitations which cannot be obtained from the orig¬ 
inal ones via the action of a local operator. We refer to 
such excitations as being in different ‘topological super¬ 
selection sectors’, or being anyons of different type; thus 
our symmetry non-trivially permutes the topological su¬ 
perselection sectors. We give a general prescription for 
understanding symmetry fractionalization in this case, 
show how it fits into the classification of SETs, and con¬ 
struct a wide class of exactly solved examples illustrating 
our results. 

One way to distinguish among different symmetry en¬ 
riched phases is based on how the symmetry fraction¬ 
alizes on the anyond^. Let us for the moment review 
the better understood case where the symmetry group G 
does not permute any of the topological superselection 
sectors. Then, because anyons are excitations that can¬ 
not be created locally, they may carry fractional symme¬ 
try quantum numbers. For example, if G is spin rotation 
symmetry SO (3), certain anyons might carry half inte¬ 
gral spins. This fractionalization of the symmetry on a 
given anyon b is captured by a collection of Berry phases 
{c Ub(g, h)}, where g, h are group elements in G: 

u^ul b) =M9,h)U^, ( 1 ) 

where U g b ^ is the ‘local’ action of g on the anyon 6, to be 
defined in more detail below. In mathematical language, 
for each anyon 6, the set of Berry’s phases {ub(g, h)} de¬ 
fines a so-called group cohomology class in iL 2 (G, U( 1)). 
The assignment of fractional symmetry quantum num¬ 
bers to anyons must also be consistent with the anyon 
fusion rules, which leads to the compatibility conditions 

u h {g, h) = u c {g, h)uj d (g , h) ( 2 ) 

whenever the anyon b is an allowed fusion product of c 
and d. 

Another way to study symmetry fractionalization is us¬ 
ing twist defects of the symmetry G® 10 ^, which are just 
extrinsic modifications of the Hamiltonian that insert a 
flux of G in a particular location. For example, an ex¬ 
trinsic defect of the Z 2 spin-flip symmetry in a nearest 
neighbor Ising model is constructed by reversing the sign 
of JcrfcrJ terms on edges ( ij ) bisected by a branch cut ex¬ 
tending from the defect core out to infinity; see appendix 


|B] for a precise definition. Such extrinsic defects have 
topological superselection sectors, which can be changed 
by fusing and splitting anyons from the defect. Suppose 
we fix a choice of such superselection sector, denoted r p , 
for each g defect. Then generically defect fusion rules 
will close only modulo an anyon ambiguity a g y. 

T g y.T h = T gh X a 9)h (3) 

Non-trivial symmetry fractionalization is then reflected 
in the fact that there is no choice of r g which makes all 
of the a 9: h = 1. These a Q: h (which we will see later can 
all be chosen to be abelian anyons) are directly related to 
the Berry phases u)b(g, h) defined above. Indeed, braiding 
the defect r g around the anyon b gives a local action U g b ^ 
of g on b discussed above, so that, using eq. [ 3 J the phase 
difference co^(g, h) in eq. [I] is just the full braiding phase 
of a 9: h around b. Because of this connection we introduce 
new notation, and denote a g ^ by uj(g,h). Thus uj(g,h), 
without a subscript, is an anyon-valued function of pairs 
of group elements, and all of the Berry phases ujb(g,h) 
can be recovered from it: 

Ub(g,h) = (u(g,h),b), (4) 

with the angular brackets denoting the full braiding 
phase. In mathematical language, w(g, h) defines a group 
cohomology class valued in the abelian any ons, i.e. an 
element of iL 2 (G, ^abelian); see appendix | b| for more de¬ 
tails. 

The discussion so far applies only to the special case 
of G acting trivially on the quasiparticle topological su¬ 
perselection sectors. To what extent does it generalize 
to a situation where the symmetry might non-trivially 
permute the topological superselection sectors? Symme¬ 
tries with such non-trivial permutation action have been 
dubbed ‘anyonic symmetriesFor example, it is pos¬ 
sible for a certain Z 2 symmetry to turn an electric ‘e’ 
exci tation into a magnetic ‘m’ excitation in the Z 2 toric 
code 13 14 . In this general permuting (or ‘twisted’) set¬ 
ting, it is difficult to make sense of fractional symmetry 
quantum numbers assigned to anyons, since even the no¬ 
tion of a local action of G on anyons does not make sense: 
e.g. there is no local operator that turns an ‘e’ anyon into 
an ‘m’ anyon in the toric code. Another complication is 
that in the permutin g case, extrinsic twist defects are 
generically non-ab elia d 8 * 11 * 1 5 -f^. This makes it more dif¬ 
ficult to study their fusion rules and extract from them 
any information about the symmetry enriched phase. 

In this paper, we study this general permuting situ¬ 
ation. Our first approach is to build concrete exactly 
solved Hamiltonians which realize symmetries that per¬ 
mute anyons. This is inspired by work of Hermele®, 
who built such models for Z n gauge theories with non¬ 
permuting symmetries, and found a class of distinct 
SET Hamiltonians naturally parametrized by a function 
l d(g, h ) in iL 2 (G, Z n ). The physical interpretation of this 
u(g,h) is exactly what was discussed above, with the 
Z n in iL 2 (G, Z n ) interpreted as the subgroup of fluxes 
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(0, fc), k = 0, 1 , ... ,n — 1 , among the set of all anyons, 
which are are just the charge j, flux fc, composites (j, fc), 
j, & = 0,1, *.., n— 1. In this special context of a Z n gauge 
theory we will abuse notation slightly and identify this 
Z n subgroup of fluxes with the multiplicative group of 
nth roots of unity e 27rlk / n , k = 0 , 1 , ..., n — 1 , so that we 
can equivalently think of uj(g,h) as being [/(l)-valued. 
This just amounts to identifying uj(g,h) with the braid¬ 
ing phase cj(i ? 0 )(g, h) of the fundamental Z n charge (1, 0) 
around u)(g,h), which contains all the information about 
u)(g, h ) in this special Z n gauge theory case. We will 
make this identification throughout sections [II] and III of 
our paper, which deal only with Z n gauge theories, and 
where it will thus not cause confusion. 


The model Hamiltonians of reference |T] are explic¬ 
itly designed to produce a Berry’s phase of u(g,h) for 
the fundamental Z n charge under the G symmetry ac¬ 
tion. Now, in our anyon-permuting situation, we find 
that we can construct a similar class of G symmetric 
Hamiltonians - again with the topological order of a Z n 
gauge theory - with only a slight modification of the 
constraints on u(g,h). These modified constraints turn 
out to define a mathematically well known generaliza¬ 
tion of group cohomology, called twisted group cohomol- 
ogy, H? wisted (G, Abelian)- The symmetry in these models 
ends up permuting the gauge charges, and also permuting 
the gauge fluxes in the same way. 

Although generalizing the models of reference [T| to the 
permuting case is rather straightforward, the physical in¬ 
terpretation of w(g, h ) is now less clear. First of all, as 
discussed above, the naive interpretation of u(g,h) in 
terms of symmetry fractionalization on the Z n charges is 
unavailable to us in this permuting setting. One can still 
study extrinsic twist defects of the symmetry however, 
and hope that u(g,h) shows up in their fusion rules, as in 
equation [3] in the non-permuting case. However, it turns 
out this is not always the case: there exist gauge inequiv¬ 
alent choices of w(g, h ) in our models which nevertheless 
give rise to the same defect fusion rules, at the level of 
superselection sectorsf 19 The corresponding Hamiltoni¬ 
ans then cannot be distinguished by the fusion rules of 
the defects, at least at the level of superselection sectors. 
Nevertheless, these Hamiltonians do define distinct SET 
phases, as we check by fully gauging G in our models and 
examining the statistics of the resulting quasiparticle ex¬ 
citations, which turn out to be different in the two cases. 
Indeed, the gauged models have the topological order of 
an E gauge theory, where E is a ‘twisted’ product of G 
and Z n , with the twist determined by uo(g,h)] distinct 
u;(< 7 , ft) give rise to distinct E. 

A more complete picture of how defect fusion data re¬ 
late to SET order can be obtained by studying defect 
fusion not only at the level of superselection sectors, but 
also at the level of ‘F-matrices’, i.e. associativity rela¬ 
tions for the fusion of defects and anyons. At this level, 
it turns out that gauge inequivalent choices of u(g, h) do 
indeed give rise to inequivalent collections of defect fu¬ 
sion and associativity data. In particular, even when the 


defect fusion rules are the same at the level of supers¬ 
election sectors for two such theories with inequivalent 
cj(g, ft), the two collections of F-matrices will be distinct 
and gauge inequivalent. In order to see this, we move 
beyond our specific class of lattice model examples, and 
develop a general framework for studying arbitrary SETs 
with permuting symmetries. The basic assumption in 
this formal algebraic approach is that extrinsic defects 
can be braided and fused with each other and with the 
anyons. Just as in the case of ordinary anyons, whose 
fusion and braiding structures - namely unitary modular 
tensor categories (UMTCs) - are highly constrained, the 
algebraic structures in the present setting involving ex¬ 
trinsic defects, so-called ‘braided G-crossed categories®^, 
are also highly constrained. Note that these are not the 
same structures, because extrinsic defects do not behave 
exactly like anyons: instead, they have branch cuts which 
are visible to the other excitations. For example, braid¬ 
ing around a defect can change anyon type, something 
that is not allowed in a UMTC. 

Classifying all braided G-crossed categories is at least 
as difficult as classifying UMTCs, since the latter are a 
subset of the former. However, in trying to distinguish 
SETs, we are really interested in the simpler problem of 
classifying all braided G-crossed categories with a given 
fixed anyon content and symmetry group G. This clas¬ 
sification problem has been solved in reference [20 and 
the resulting mathematical machinery has been applied 
to classify SETs in references [ 8 ] and 0 Using this general 
classification, one finds an invariant which distinguishes 
braided G-crossed categories with the same permutation 
action of G which is valued in H 2 wisted (G , M a beiian)- This 
invariant reduces to the ordinary fractionalization class 
in the non-permuting case, where it is seen in the de¬ 
fect fusion rules already at the level of superselection sec¬ 
tors. In the more general permuting case, though, it can 
generically only be obtained from knowledge of both fu¬ 
sion rules and F-matrices involving 2 defects. We will 
review the stepwise construction of braided G-crossed 
categories, following reference [ 20 ] and using an intuitive 
graphical calculus, and see explicitly how the invariant 
shows up in fusion and F-matrix data. 

To connect this formal approach to our class of model 
Hamiltonians, we study a specific example: a Z4 gauge 
theory with a Z 2 symmetry acting by (z, j) (4 — i, 4— j) 
on the charge/flux composites. There are two distinct 
lattice Hamiltonians of the type we consider with this 
symmetry action, corresponding to two inequivalent sets 
of Berry phases u; + (g, ft) and cj_(g, ft), and they are ex¬ 
actly of the type discussed above: their defect fusion rules 
are identical at the level of superselection sectors, but 
they correspond to distinct SETs, because they gauge to 
different topologically ordered theories. Therefore, they 
should differ in their F-matrix data, for F-matrices in¬ 
volving two defects and an anyon. It is difficult to ex¬ 
tract such F-matrix data from the lattice Hamiltonians 
directly, but fortunately, because of the strong algebraic 
constraints within the braided G-crossed category, this 
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F-matrix data is also reflected in defect braiding data. 
Specifically, we will see that, for this example, the F- 
matrix data should be encoded in certain anyon-defect 
braiding processes, and we explicitly confirm that this is 
the case for our lattice models. 

The remainder of the paper is structured as follows. 
In section [IT] we construct our exactly solved lattice SET 
models. As in reference[lj they are given by coupling |G| 
copies of a Z n gauge theory (here \G\ is the number of el¬ 
ements in the group G ), although in our case the symme¬ 
try action non-trivially permutes the Z n -charges among 
themselves, and similarly for the Z n -fluxes. In section [III 
we explicitly gauge G in these models, study the topolog¬ 
ical superselection sectors of their defects, and show that 
the topological order of the gauged theory is the quantum 
double of the non-central extension of G by Z n deter¬ 
mined by w(g,h), generalizing the non-permuting result 
of 1 . In particular, whenever these non-abelian gauge the¬ 
ories are distinct, so are the underlying SETs, showing 
that this construction does indeed produce non-trivial 
SETs. Of course, these SETs are far from the most gen¬ 
eral ones possible - in particular, since after gauging we 
obtain discrete gauge theories, all of our defects have in¬ 
tegral quantum dimension. Nevertheless, they still form 
a wide class of explicit realizations of the various phases 
allowed by the recent general classification of SETs in 
two dimensions. In particular, we discuss in detail the 
simplest example of Z 4 gauge theories with symmetry 
G = Z 2 acting by k — k for k G Z 4 , where there are 
two symmetry enriched phases, which give non-abelian 
B 4 (dihedral group of symmetries of the square) and 
Qs (quaternion group) gauge theories respectively upon 
gauging G. Finally, in section IV we develop the general 
theory of defect fusion rules and their deformations, ap¬ 
plicable both in the non-permuting and permuting cases. 
We use a graphical formalism to introduce the mathe¬ 
matical description of defect superselection sectors, and 
describe defect fusion rules within this formalism. Math¬ 
ematical results of reference SOI then allow us to enumer¬ 
ate the gauge equivalence classes of such defect products, 
and show that they are in one to one correspondence with 
^twisted (G, A). We then again study the Z 4 gauge theory 
example mentioned above, and treat it within the con¬ 
text of this general theory. Finally, we summarize and 
discuss new directions in section El 


We will take G to be abelian for convenience, though 
we believe our results generalize to non-abelian G. Al¬ 
though we work with the topological order of an abelian 
Z n gauge theory, our results readily generalize to ar¬ 
bitrary abelian groups. We also treat the special case 
G = Z 2 , n = 4 in detail. 


A. Hilbert space 

Our model is a Z n gauge theory living on a certain 
oriented, quasi-2d lattice. Following reference[lj we start 
with a truly 2 d oriented lattice, which can be taken to be 
a square lattice in the xy plane for all of the examples we 
consider, and stack \G\ identical copies of it. This stack¬ 
ing allows us to identify corresponding vertices and links 
in each copy. In particular, consider the set of \G\ vertices 
that all have the same x, y coordinate. For any ordered 
pair (i?, w) of such vertices, we add an oriented link con¬ 
necting v to w. For clarity, we refer to these |G|(|G| — 1 ) 
links as vertical , as opposed to the links within layers, 
which will be called horizontal The orientation of hori¬ 
zontal links is the same across all layers. 

The set of |G| corresponding vertices together with the 
|G|(|G| — 1 ) links connecting them will also be referred to 
as a supervertex (reference (Tj calls this a Cayley graph). 
Likewise, the set of |G| links which project to the same 
2 d link will be referred to as a superlink. 



FIG. 1: Lattice on which our model is defined. The 
links carry Z n labels. 


II. EXACTLY SOLVED LATTICE 
HAMILTONIAN 

In this section we write down a family of exactly solved 
lattice models of G-symmetric Hamiltonians, with G per¬ 
muting the anyons. The goal here is simply to describe 
the Hilbert space, operators, and symmetry action in as 
explicit a way as possible, and motivate the form of the 
Hamiltonian in equation [37] In later sections we analyze 
the models described by this Hamiltonian in detail, and 
see that they correspond to distinct SETs. 


The Hilbert space is taken to be spanned by Z n la¬ 
bellings of the links of our lattice. From now on we will 
identify Z n with the set of n’th roots of unity, i.e. com¬ 
plex numbers of the form e 27Tl j/ n ^j = 0,..., n — 1. For¬ 
mally, we define an n dimensional link Hilbert space Hi 
whose basis states are in one to one correspondence with 
such roots of unity: 

Hi = span{ | rji) \rji G {e 2mj/n ,j e Z}} (5) 

and take the total Hilbert space H to be the tensor prod¬ 
uct of these link Hilbert spaces (including both horizontal 
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and vertical links): 


example of a G -invariant Hamiltonian is: 


U = (g)Ui. (6) 

l 

On each link Hilbert space we define the usual ‘phase’ 
and ‘charge’ measuring operators ai and ep. 

ei\m) = \e 2 * i/n m) ( 7 ) 

ai\m) = mh) ( 8 ) 

By tensoring with the identity on all other links, we can 
think of ai and ei as being defined on the total Hilbert 
space 1~L. Note that two such operators acting on different 
links Z and V commute. 

The Hamiltonians we will work with contain terms 
which act on certain groupings of links, associated to ver¬ 
tices and plaquettes, and before we can write them down 
we need to establish some effective notation. First of all, 
as we mentioned above, our quasi- 2 d lattice is oriented, 
which means that there is a preferred choice of direction 
for each link. This orientation is efficiently encoded in a 
function s v (l ), where Z is any link and v is one of the two 
endpoint vertices of this link: 


s v (Z) 


1 if Z points toward v 
— 1 if Z points away from v. 


( 9 ) 


We will assume that our orientation is consistent across 
the \G\ layers, i.e. s v (l) = s v /(l') whenever v,v' are in 
the same supervertex, and horizontal links Z, V are in the 
same superlink. 

Additionally, we now also assign an orientation to all 
plaquettes p (i.e. plaquettes involving any combination 
of horizontal and vertical links). This orientation is just 
a choice of direction, either clockwise or counterclock¬ 
wise, along the links that border p. For each such link Z 
bordering a plaquette p, this choice of direction could be 
the same or opposite to the one defined by Eq. [9j This 
distinction is encoded in a function s p (Z), where Z is a link 
bordering the plaquette p: 


s p (0 _ 



if Z oriented with p 
if Z oriented against p 


( 10 ) 


We will assign this plaquette orientation consistently 
across the layers, in that if p and p' are plaquettes made 
up entirely of horizontal links that project to the same 
plaquette in the xy plane, we assign them the same orien¬ 
tation. This just means that if Z and l' are corresponding 
links of p and p' respectively (so that Z, l' are in the same 
superlink), then s p (l) = s p '(l'). Besides this constraint, 
the plaquette orientations are chosen arbitrarily. 

Now, in reference [Q G acts by permuting layers, and 
since such a permutation induces a one to one mapping 
of the underlying oriented quasi- 2 d lattice to itself, an 


w 

1 

II 

S3 

- Bp + h.c. 

( 11 ) 

V 

Av= n e ? o(,) 

l~v 

b p = IK pa) - 

lep 

p 

( 12 ) 

(13) 


where the notation Z ^ v refers to all links Z that begin 
or end at v, and Z G p refers to all the links that border a 
given plaquette p. While eq. [Tl] is adequate in the case 
where the symmetry does not permute the gauge theory 
quasiparticles (anyons), we will need a slightly different 
construction for a symmetry action which does permute 
the anyons. 


B. Symmetry action and Hamiltonian 

In our model, G will act both by permuting the links 
and changing the Z n labels on the links. The permuta¬ 
tion of the links induced by g G G is the same as that in 
reference [1: given a vertex v in layer h, we let gv denote 
the vertex in layer gh which is in the same supervertex 
as v. Then, for a link Z = (vv'), we define gl = ( gv gv'). 
The change in the Z n label that goes together with this 
link permutation - which is the new feature of our model, 
and is referred to as a twisting - is encoded in an inte¬ 
ger valued function p(g), with each p(g) relatively prime 
to n (i.e. having no common factors with n). p is re¬ 
quired to satisfy p(gh) = p(g)p(h ) (note that this is mul¬ 
tiplication of integers modulo n) and allows us to define 
a permutation action of G on Z n , namely g g p( < 9 \ 
More explicitly, if g = e 27rz/c / n , then this action just takes 
k —>• p(g)k mod n. An example that we will focus on in 
the rest of the paper is G = Z 2 and n = 4, and p(g) = — 1 
for the non-trivial generator g of Z 2 . 

Using p, we define the global action of G on the Hilbert 
space as follows. With a slight abuse of notation, we 
denote the unitary action of g by U g , regardless of what 
Hilbert space is being acting upon. For the link degrees 
of freedom we let: 


II 

to 

(14) 

This induces the action on operators: 


© 

II 

© 

ccT 

1 

(15) 

U g eiUp=e p g \ a) 

(16) 

We can immediately infer that the action on 
plaquette terms defined in eq. [Tl] is 

vertex and 

U g B v Up = B^- 1 ) 

(17) 

U g A v Up = A$\ 

(18) 
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where gp is the plaquette made up of the links gl , for all 
l G p. Note that for p(g) ^ 1, the Hamiltonian defined 
in eq. 11 is not invariant under this global action of g. 
Instead, we consider the following more general Hamilto¬ 
nian: 


To construct the plaquette terms, we take a trace of the 
matrix produced by a closed loop of a/’s. For example, 
consider a triangular plaquette p with links 1,2,3, and 
uj p = 1 , and suppose |T) is an eigenvalue rjj eigenvector 
of each ctj. Then: 


h set = - E (e^™ + E• (is) 

m= 1 \ v p / 

Here uj p are phases - in fact, nth roots of unity - associ¬ 
ated with each plaquette p, which satisfy 


fjj — / .PCs 1 ) 

LUgp — UJp 


( 20 ) 


We can verify that Hset is G-invariant: 


UgUsETUp = 


n / 

- E E- 4 ™ 0 ”+E 

772—1 \ V p 

n / 

= - E E- 4 ™ 0 ”+E 


772=1 


= H SET, if U) gp = U>M 


( 21 ) 


19 


describes a 


Thus, for l dp which satisfy u gp = ujp" 9 \ eq. 

Hamiltonian that is invariant under the twisted G action. 

Let |T) be a state of B corresponding to a specific 
Z n labeling of links. Recalling that the spectrum of B p 
consists of the roots of unity e 27Tlk / n , we see that 


EK-%)”» 

772=1 


n|M/) ifB p |tt)=u; p |tt> 
0 otherwise 


Thus the operator defined on the left side of eq. [22] is 
equal to n times a projector. We now describe a no¬ 
tation that will let us concisely express such operators; 
we emphasize that this formulation is nothing more than 
a notational convenience. First, recall that the regular 
representation of a group H is an \H\ dimensional vector 
space with basis {\h')\h f G H}, where h G H acts by 


h : | ti) | hti) (23) 


Thus h G H is represented by an \H\ by \H\ matrix M(/i), 
where each column and row is labeled by a group element 
g,k G iL, and the matrix elements are 


M(h)g tk = Sg^hk (24) 

A feature of these matrices is that Tr M(h) = \H\5h,i> 
Now, recall that the operator a/ acts by the phase pi on 
1 77 /). In our new notation, acting with the operator a/ 
yields the matrix M(rj), where now H = Z n . In the case 
where n = 4 we have 


M 



/0 0 0 1\ 
1 0 0 01 
0 10 0 
Vo 0 1 0/ 


(25) 


(Tr ai a 2 a 3 ) |$) = (Tr MfaJMfoJMfas)) |^) (26) 



if V 1 V 2 V 3 = 1 
otherwise 


(27) 


Thus, the operator in 22 can be rewritten in this new 
notation as: 


E ^pB p ) m -> Tr copBp (28) 

772=1 

yielding a notationally convenient way of writing n times 
the projector onto the eigenvalue 1 subspace of uj~ x B p . 
Here it is understood that the complex number uj p is 
substituted with its regular representation matrix. No¬ 
tice that the trace on the right hand side of this equation 
is over the auxiliary regular representation, and not the 
many body Hilbert space; both sides are operators in the 
many body Hilbert space. 

With this new notation, our Hamiltonian takes the 
form 

H set = -E(E^)-E TV < B p ( 29 ) 

One benefit of our new notation is that it makes it easy 
to generalize to the case of an arbitrary abelian 

gauge group iL, rather than just Z n . Indeed, to do this 
one just needs to let p be a map from G to Aut(iJ), the 
group of automorphisms of H. Nevertheless, we will stick 
to a Z n gauge group in the remainder of this paper. 


C. Supervertices and superlinks 

The next step is to discuss the phases cj p , different 
choices of which will give rise to different twisted sym¬ 
metry enriched phases. To facilitate this discussion, it is 
useful to first group the degrees of freedom in our model 
in a slightly more convenient way. First, recall that a 
supervertex V is a collection of \G\ vertices which all 
project onto the same point in the xy plane, i.e. are ver¬ 
tically aligned. We now tensor the Hilbert spaces living 
on these links into a single supervertex Hilbert space By 
spanned by {\py (g, ft))}, where r]y(g,h) is the Z n label 
of the link l g , g h connecting layer g to layer gh (here h 
is necessarily different from the identity in G ). We then 
denote by ay(g,h),ey(g,h) the action of the operators 
a i g gh 5 e i g gh 5 tensored with the identity on the remaining 
links, in the Hilbert space By. 

Similarly, we define a superlink L to be the collection 
of \G\ horizontal links whose projections to the xy plane 
are all the same, and likewise define a superlink Hilbert 
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space 1~Ll to be the tensor product of the associated |G| 
link Hilbert spaces. We denote by a^g), ez,(g) the action 
of the operators ai g ,ei g on the link l g G L in layer g, 
tensored with the identity on the remaining \G\ — 1 links. 

Our total Hilbert space H is thus a tensor product of 
the supervertex and superlink Hilbert spaces: 


n= ((g) ® ((g) 


(30) 


D. Conditions on the fluxes uj p 

We now discuss the choice of U(1 ) phases c jJ p , which we 
also refer to as Z n fluxes, since they are restricted to take 
values in the n’th roots of unity. First of all, through¬ 
out this paper we will deal exclusively with the situation 
where the only non-trivial uo p (i.e. uj p ^ 1) occur for pla- 
quettes p that sit entirely within a single supervertex, or, 
in other words, contain no horizontal links. Let us focus 
on a specific supervertex V. Each vertical plaquette p 
within it is labeled by a triple (/, g,h), where g,h ^ 1, 
indicating that it involves the links connecting layers /, 
fg and fgh. The corresponding term in the Hamiltonian 
(eq. [29| is: 


Tr u p 1 B P M Tr uj p 1 B v (f,g , h ) 

= Tr wpavif, gh)~ 1 a v {fg , h)a v (f, g) (31) 

While most plaquettes will be 3-edged, we are allowed to 
set h = g -1 , producing a degenerate, 2-edged plaquette. 
This situation can still be captured by eq. [31] by defining 
%(/, e) = I (the identity operator). 

Using eq. [20] we see that all of the uj p are uniquely 
determined by the uj p for p = (e,g, h), i.e. for plaquettes 
p which start at the identity element of G. Letting u ;(#, h) 
denote uj p for p = (e, g , h), we then have that for a general 
plaquette p' = ( f,g , h), ay = w(g, h) p G\ 

f9 



FIG. 2: A diagram showing the flux piercing the 
plaquette (/, g, h). Note that it spans layers /, fg and 

fgh 


Now, we would like to work with models which are 
unfrustrated, i.e. whose ground states are lowest energy 
eigenstates of all of the plaquette terms. Consider the 


f -a fg 



FIG. 3: A tetrahedron spanning the layers 
{e, /, fg, fgh}. The total 7L n flux emanating out of this 
tetrahedron must be trivial (eq. [32| in order to avoid a 
degenerate set of frustrated ground states. 


tetrahedron formed by the layers {e, /, fg, fgh}, which 
contains 4 plaquettes: (e,f,g), (f,g,h), (e,fg,h), and 
(e,f,gh). A necessary and sufficient condition for the 
model to be unfrustrated is that the 7L n flux emanating 
out of any such tetrahedron be zero, i.e. 

u(g, h) p ^uj{f, gh) = w(/, g)u(fg, h). (32) 

Indeed, if the ground state is unfrustrated, there must 
be some labeling {py (g, h)} of the links in the superver¬ 
tex V (namely one that corresponds to a configuration 
that enters the unfrustrated ground state with non-zero 
amplitude) such that 

'nv{f-,g)'nv{fg,h)riy 1 {f,gh) = ui(g,h) p{f) ,Wf,g,h e G 

(33) 


Expressing uj in terms of the rjy using this equation, we 
see that it satisfies eq. JH Conversely, given a choice 
of cTs that satisfy eq. 32, we can simply set gy(g,h) = 
w(g,h). This link labeling satisfies all of the plaquette 
terms within the supervertex V. Later, we will see that 
it can be extended to an unfrustrated ground state of all 
of the vertex and plaquette terms in our model - indeed, 
we will explicitly solve the model for any choice of fluxes 


satisfying eq. 32 


Certain different choices of c o(f,g) actually define 
Hamiltonians which can be made equivalent by redefining 
link variables: 


Vv{g, h) ->■ p(h) p{9 \ v {g, h) (34) 

where g : G —* Z n is an arbitrary function. This redefi¬ 
nition then takes 

w(ff, h) u'(g, h) = u>{g , h)g(h) p(a) g{g)g{gh)~ l . (35) 

Note that the new c o'(g,h) also satisfy equation [32] In 
group cohomology terms this means that equivalence 
classes of non-frustrated Hamiltonians of the above form 
are parametrized by twisted group cohomology classes in 
H 2 p (G, Z„)P 

For simple enough G, it is easy to compute these co¬ 
homology groups explicitly. For example, take G = Z 2 = 
{ 1 ,- 1 }, n = 4. Then there is only a single plaquette, 







bounded by the links 1 -1,-1 —^ 1, which is pierced 

by a Z 4 flux 1 , —1), subject to a twisting p(—1 ) = — 1 . 
Equation [32] produces a non-trivial constraint only for 
f,9,h= -1: 

w(-l,-l)^>=w(-l,-l) (36) 

which implies u(— 1 , — 1 ) = ± 1 . As we will see, these 
two choices of a;, which we call lj±, will produce two 
inequivalent SET Hamiltonians, which in turn yield two 
distinct non-abelian gauge theories once we gauge the Z 2 
symmetry. 

E. Form of the SET Hamiltonian 


Note that there are two distinct kinds of vertical plaque- 
ttes: ones contained entirely in a single supervertex, and 
ones involving a superlink and the adjoining two super¬ 
vertices. Only for the ones contained entirely in a single 
supervertex can we have uj p 7 ^ 1 . 


v 2 p v 2 



Let us write out the final form of the SET Hamilto¬ 
nian in a form convenient for gauging G, which we do in 
the next section. Recall (eq. 30) that our total Hilbert 
space T-L is a tensor product of supervertex and superlink 
Hilbert spaces Hy and T-Lp. Because all of the links in 
a superlink are oriented the same way we can define the 
orientation factor sy(L) = s v (l) = ±1 where u, l are any 
vertex and adjoining link in V and L respectively. Now, 
plaquettes containing only horizontal links can similarly 
grouped into superplaquettes P. Again, because the pla- 
quette orientations have been chosen so that all plaque¬ 
ttes p in a superplaquette P are oriented the same way, 
we can define sp(L ) = s p (l) for any p in P and l in L 
bordering p. We then have the form of the Hamiltonian: 


FIG. 4: Vertical links are colored red and horizontal 
links are colored black. Collections of vertical links 
connecting points with the same x-y coordinate 
constitute supervertices, labeled V\ and W Plaquettes 
made entirely out of horizontal links that project to the 
same plaquette form a superplaquette, labeled P in the 
figure. Plaquettes containing at least one vertical link 
are deemed vertical. The dotted lines identify vertices 
and help offset the figure. Certain links are dashed in 
order to add perspective. 


III. DISTINGUISHING SET PHASES 


iw = -e (e( E ^foA + E 1 * 5 ^)) 

geG \ V \mez n / p J 

- E Trw-'Ci, (37) 

vertical p 

where V and P range over supervertices and superpla¬ 
quettes respectively, and the various terms in the above 
sum are defined as follows. We’ve rewritten our vertex 
and plaquette operators as 

Av{g) = H e v (g,h)~ 1 e v (gh,h~ 1 ) 

heG 

X n e L(g) Sv(L \ (38) 

L~y 

which denotes the term acting on the vertex on layer g 
of the supervertex V , and 

Bp(g)= n a L {L) (g) (39) 

lep 

denotes the term acting on the plaquette on layer g of 
the superplaquette P. Finally, for any plaquette p which 
is not composed solely of horizontal links, which we refer 
to as ‘vertical’ in eq. [37] above, 

q> = IK p(,) - ( 4 °) 

lep 


In the case of symmetries that do not permute anyons, 
the fluxes w(g,h) appearing in the Hamiltonian in equa¬ 
tion [37] can be physically interpreted as Berry phases 
for the symmetry action on the fundamental Z n charg^. 
However, such an interpretation does not generalize read¬ 
ily to the anyon permuting case, and this motivates us 
to couple the model to a G gauge field and examine the 
resulting ‘gauged’ theory. There are two different ver¬ 
sions of such a gauged theory: one can either make the 
G gauge field a fully dynamical degree of freedom, or one 
can treat it as a background probe field. 

In the case of a dynamical gauge field, it turns out 
that different SET Hamiltonians - i.e. different choices 
of h ) - can be distinguished by the statistics of the 

excitations in the gauged theory. Demonstrating this fact 
will take up the bulk of this section. Indeed, to fully un¬ 
derstand the gauged theory with dynamical gauge field 
G, we first perform a minimal coupling procedure of the 
Hamiltonian in eq. [37] to such a dynamical G gauge field, 
and then perform a series of transformations, analogous 
to those in reference \l\ to simplify the form of the result¬ 
ing gauged Hamiltonian, without altering the low energy 
physics. Although these transformations are technically 
complicated, there is a simple intuitive picture for what 
is going on: essentially, the gauge field G inserted along 
superlinks should be viewed as allowing permutations be¬ 
tween the different layers. When the symmetry is gauged, 
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any potential physical distinction between the different 
layers is therefore lost, and hence the physical states in 
the gauge theory live on an ordinary 2 d lattice, as op¬ 
posed to a |G|-fold stacked quasi- 2 d lattice. Indeed, we 
find explicitly that the simplified gauged theory is just 
a discrete gauge theory of a group E on an ordinary 2d 
lattice. Here E is a group which has Z n as a normal sub¬ 
group, E/ Z n = G, and p(g)rj = g~ l gg , where g is any 
lift of g in E. E is called the group extension of G by Z n 
determined by the permutation p and an element of [a;] 
of i?2(G,Z n ). Each uj corresponds precisely to one such 
group extension (See Appendix [d]). 

The second version of a G gauged theory is one where 
the G gauge field is treated as a background field - we will 
refer to this as the non-dynamical case. Here the fluxes 
of G are not dynamical excitations, but rather extrinsic 
defects, requiring a branch cut in the Hamiltonian. One 
reason one might want to examine this case is that it 
is conceptually simpler than that of the fully dynamical 
G gauge field, as it requires no extra degrees of free¬ 
dom. Another reason is that recent work (reference [27]) 
has studied the general mathematical structure encoded 
in such extrinsic defects, called a braided G-crossed cate¬ 
gory. We will also discuss general aspects of such braided 
G-crossed categories later in the paper, but in this sec¬ 
tion we will just study them in the context of the class of 
lattice models we have just introduced. Having already 
coupled these models to a dynamical G gauge field, it 
turns out that the analysis of the non-dynamical case is 
easy: in the final form of our dynamical gauged model as 
an E gauge theory on a 2d lattice, it just amounts to in¬ 
cluding only vertex terms corresponding to Z n C F 1 , and 
setting the coefficients of the vertex terms corresponding 
to other elements of E to 0. 


One may wonder how two of our SET models corre¬ 
sponding to two distinct choices of [cj] G iT^(G, Z n ) can 
be distinguished when coupled only to a non-dynamical 
G gauge field. In particular, can different choices of [u] 
lead to different extrinsic defect types? It turns out that 
the answer is no: the properties of a single extrinsic g- 
defect are uniquely determined by g and the permuta¬ 
tion p, and independent of [cj]. Specifically, the set of 
topological superselection sectors bound to a single ex¬ 
trinsic p-defect, the fusion rules of these superselection 
sectors with external anyons (that is Z n gauge charges 
and fluxes), and the F-matrix associativity constraints 
involving two anyons and a defect are all uniquely deter¬ 
mined by g and p (the F-matrices are unique only up to 
the appropriate gauge degree of freedom). This collec¬ 
tion of data is known as an invertible bimodule category, 
and will be discussed in section [IV[ Thus, in order to dis¬ 
tinguish SETs with different [a;], we have to go beyond 
the case of a single extrinsic defect, and look at pairs 
of such extrinsic defects and their fusion rules. Indeed, 
in the non-anyon-permuting case of reference |l] different 
choices of [a;] lead to different fusion rules for extrinsic g 
defects, at the level of topological superselection sectors. 

In our anyon-permuting case, the situation is more sub¬ 


tle. It is still true that different choices of [uj] lead to 
different fusion rules for extrinsic g defects, but not at 
the level of topological superselection sectors. In other 
words, the fusion rules (which tell us which topological 
superselection sectors in defect gh can end up as the fu¬ 
sion product of specific sectors in defects g and h) might 
be the same for different choices of [cj]. In this case, the 
distinction between two such different choices of [c u\ can 
only be seen in the F-matrices involving two defects and 
an anyon. More precisely, since the F-matrices are not 
gauge invariant, the distinction can only be seen in the 
gauge equivalence classes of such F-matrices. We explain 
this precisely in section [IV C| below. 

We will discuss the general theory of defect fusion in 
section |IV[ At the end of the present section, however, 
we will analyze the specific example of a Z 4 gauge the¬ 
ory with G — Z 2 = {1, —1}, with p{— 1 ) = — 1 . We will 
see that iT^(G, Z4) = Z2, so there are two inequivalent 
choices of u; (/,#), which we denote c and c j_. The 
group extension E corresponding to the trivial co-cycle 
is Bg (the dihedral group on 4 points), while that 
corresponding to the non-trivial one cj_ is the quater¬ 
nion group Qg- We will find that, regardless of whether 
we choose c or cj_, the extrinsic Z 2 defects in the two 
SETs have the same superselection sectors and the same 
fusion rules for these sectors: as discussed above, these 
are independent of u). Thus, to tell the difference between 
the two theories, we must probe more subtle data. In¬ 
deed, we can either look at the quasiparticle statistics in 
the dynamical G gauged theory, i.e. the Bg and Qg gauge 
theories, and see that they are different, or, as discussed 
above, we can detect the difference in the F-matrices cor¬ 
responding to the defect fusion rules. 


A. Gauging prescription 

The goal of this somewhat technical appendix is to de¬ 
rive e qs. [5Q] and[57l which describe our SET Hamiltonian 
in eq. |37| coupled to a G gauge field. First we introduce 
G gauge field degrees of freedom, which are just |G| di¬ 
mensional Hilbert spaces which we insert between any 
supervertex and an outgoing superlink (see figure [TT]) . 
The |G| degrees of freedom are incorporated by enlarg¬ 
ing each superlink Hilbert space El —> El ® C^l, with 
states in this larger Hilbert space carrying an extra G 
gauge field label g^ 

\{vL(g)}g) -> \{v L (g)}g,gL) (41) 

The minimal coupling prescription we use is as follows. 
Given any local term in the original ungauged Hamilto¬ 
nian, for a G gauge field configuration which is gauge 
equivalent to the trivial configuration in the vicinity of 
this local term (i.e. has no G fluxes), the form of the cor¬ 
responding minimally coupled term is completely fixed 
by G gauge invariance. For a G gauge field configura¬ 
tion which does contain nonzero G fluxes in the vicinity 
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we first define an auxiliary ‘book-keeping’ operator 


aL(\{nL(g)}g,9L)®---) = 

m ( 9 l)U ~ l 1 {\{rt L {g)}g,g L ) ® • • •), (44) 


where M{gi) is the \G\ by \G\ matrix representing in 
the regular representation of G (see the discussion around 
eq. 23), U Ql is the global action of gL (as defined in eq. 


14), and the ellipses denote the rest of the Hilbert space 


P. Now define 


FIG. 5: Placement of the G gauge variable relative to 
the supervertex and superlink 


Bp(g) = n K(sK) sp(i) • (45) 

lep 


of this local term, we simply set the corresponding min¬ 
imally coupled term to 0. This actually only occurs for 
the superplaquette term, when there is a non-zero G flux 
through it. Since the Hamiltonian consists of commuting 
terms which all have negative eigenvalue on the ground 
state, setting this minimally coupled term to 0 is actually 
an energetic penalty for the non-zero G flux. Lastly, we 
include ‘vertex’ terms which make the G gauge field fluc¬ 
tuate and thus energetically impose G gauge invariance. 
These vertex terms are the only ones which alter the G 
gauge field configuration. 

The above paragraph specifies the minimally coupled 
Hamiltonian uniquely, but to actually write it out in a 
compact form it is useful to introduce some additional no¬ 
tation. It is easiest to start with the vertex terms, given 


by eq. [38] A particular such term involves a specific ver¬ 
tex v, which is part of a supervertex V. To minimally 
couple it, we have to modify each e^g) term correspond¬ 
ing to an outgoing superlink L from V (that is, one with 
sy{L ) = —1) to take account of the gauge field gL, and 
replace it with: 


e L {g) R \{vL{.g')} g’, 9 l) = 


I {VL(g')e 


— 27rz 


p(g l) 


' 9 ’ 9 ' 9L }g',9L ) (42) 


Here R is just a superscript. To keep the notation com¬ 
pact, we define a superscript-valued function s' v (L ) = R 
when sy{L) = — 1, and s v (L) being trivial otherwise, so 
that we can write the minimally coupled version of Ay{g) 
simply as: 


Av{g) = JJ e v(g,h) 1 ey(gh 1 h x ) 
h£G 

x n e L(g) s ' v(L) (43) 

L~V 

Now let us minimally couple the horizontal plaquette 
terms Bp(g). When there is a non-trivial G flux through 
P, we simply set Bp(g) = 0. When this flux is trivial, 


Note that even though the operators oll are non-local 
(because they involve the global U 9l operator, when the 
G flux through P is trivial their combined action in 
Bp(g ) cancels away from the plaquette P, and we end 
up with a local operator. Also recall that it is the trace 
of Bp(g) over the auxiliary regular representation space 
which appears in the Hamiltonian. 

Next, let us minimally couple the term C p which in¬ 
volves both a superlink and the adjacent two superver¬ 
tices. Any such p is a rectangle, and C v is a product of 
4 link terms: 


Cp = (a L (gh)a v (g,h)) 1 a v >(g,h)a L (g) (46) 


Minimally coupling this inserts the gauge field, in the 
form of the oll operators, along the horizontal links: 


Cp = (a L (gh)a L a v (g,h)) 1 a V '(g,h)a L (g)a L (47) 


Recall that there is never any G flux through such a pla¬ 
quette p. Finally, the minimal coupling of C p for a pla¬ 
quette p entirely within a supervertex V does not involve 
the gauge field at all, C v = C p . 

In summary, the minimal coupling is done by starting 
with the Hamiltonian in eq. [37] and making the modifi¬ 
cations 


az,(s) -s- a L (g)a L 

if L V 

(48) 

sv{L) ->■ Sy(L) = 

if L <- V, 

(49) 
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resulting in a minimally coupled Hamiltonian H m . c .: 


H n 


E £ E A v(s) +e^^( 9 ) 


geo 




- E Tr d> 

vertical p 

Av(g) = Yle v (g,h)~ 1 e v (gh,h~ 1 ) 
heG 

x n e L{9y' v{L) 

L~V 

Bp(g) = n Mg)a L y p{L) 

LeP 

'uj- 1 a v (f,gh)- 1 a v (fg,h)a v (f,g) 

^ ^ U if p is in V 

P | ( a L {gh)a L a v {g,h))- 1 a V ’(g,h)a L (g)a L 
K if p is in L = (' VV') 


(50) 


(51) 

(52) 


(53) 


Here the traces are over both G and Z n . The Hamilto¬ 
nian H m . c . describes our SET in a fixed background G 
gauge field configuration, corresponding to some set of 
extrinsic defects. We will also want to have a Hamilto¬ 
nian where the G gauge field is dynamical. To obtain 
it, first define operators which change the value of on 
a single superlink L. We will need two such operators, 
e L (g) and e%(g): 

eLmvL(g')} g ’,gL) = | {vLig-'g'y^h^ggL) (54) 
eL(g) R \{vL(g')}g',gL) = I (55) 

Using these we define: 


■Av(g) = u~\v) n £dgY' v(L) ) (56) 

where U g (V) is the action of g on the supervertex V. The 
Hamiltonian with dynamical G gauge field then becomes 

Hg aU g e d = H m . c . ^ ^ ^ ^ A v (g) (57) 

V geG 

This is the model whose topological order we have to 
analyze. 


B. Analysis of the topological order 


We claim that the gauged model defined by eq. 57 


is 


equivalent to - i.e. has the same topological order as - an 
E gauge theory, where the finite group E is a particular 
extension of G by Z n . We derive this equivalence care¬ 
fully in appendix [Aj while in this section we just write 
down the result. First, let us discuss the extension E. 
As a set, it is just the product Z n x G, so we can label 
its elements by ( 77 , 0 ), with g an n’th root of unity and 


geG. On the other hand, the group multiplication law 
is determined by both by p and uj(g,h): 


(vi,gi) • (772,02 ) = (t?i(772) p(si) w(5i,5 2 ),0i02) (58) 

Now we write down an E gauge theory on a single copy 
of the 2 d lattice we considered above (as opposed to 
\G\ stacked copies of it). Since the vertices and links 
of this 2 d lattice are in one to one correspondence with 
the supervertices and superlinks of the \G\ stacked lat¬ 
tice, we will just label them by V and L respectively, and 
keep calling them supervertices and superlinks to avoid 
confusion. Then the degrees of freedom are just pairs 
xl — (with 7]l an n’th root of unity and gL G G, 

and x l treated as an element of the group E) defined on 
the superlinks L. The Hamiltonian consists of two kinds 
of terms, vertex terms and plaquette terms. 

Intuitively, we define the vertex term at V correspond¬ 
ing to x G E by multiplying all of the link variables 
on links L terminating at k by r. However, there is a 
complication due to the fact that some of the links are 
outgoing and some are incoming: in one case we want 
to left multiply by x and in the other we want to right 
multiply by x ~ x . We handle both cases by defining the 
link multiplication operator: 


e v ,L(x)\x L ) = \xx L ) if sy(L) = 1 (59) 

= \x L x ~ 1 ) if s v {L) = -1 (60) 

With it we then define 

Ay(x) = ( e V,L(x )) (61) 

L~V 

The plaquette term Bp associated to any (su¬ 
per) plaquette P is defined by assigning an energetic 
penalty for E flux through P. More formally, in the 
basis of link labellings, we define: 


Bp = 1 if JJ xp Sp ^ = e 
LeP 

= 0 otherwise. 


(62) 

(63) 


Using these operators, the Hamiltonian becomes: 

h e = -EE^)-E^ ( 64 ) 


v xeE 


The equivalence between the Hamiltonian in eq. 57 and 
that in eq. [64] is non-trivial, and derived in appendix 
[A] Here we will just make a couple of comments about 
how these two Hamiltonians are related. First of all, the 
G-flux through any plaquette, given by 


n 3 l sp(l) 

lep 


( 65 ) 
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in the notation of eq. [57] is given by the projection from 
E to G of 


n 

LeP 


XL 


sp(L) 


( 66 ) 


If we represent each xl as xl = (pzo#l), then this just 
reduces to the expression in eq. [65] 

Now let us examine configurations with trivial G gauge 
field, i.e. ql — 1 f° r all L. Then the Hamiltonian in eq. 
[57] just reduces to the original Z n gauge theory on the 
\G\ fold stacked lattice. On the other hand, the E gauge 
theory given in eq. [64] reduces to a Z n gauge theory on 
the ordinary 2d lattice. The identification between these 
two is non-trivial, but in particular a Z n flux in the latter 
corresponds to that same flux penetrating all \G\ layers 
in the stacked lattice. 

Now we will examine some consequences of this equiv¬ 
alence, for a particular example. 


C. Example: Z 4 gauge theory 
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FIG. 6 : Z 2 Gauge field configuration with two widely 
separated Z 2 defects. The Z 2 gauge field is non-trivial 
on the vertical links intersecting with the dashed blue 
line. We imagine acting with e/ (the operator which 
increments Z 4 labels) to the power labeled in the figure 
(1 or —1) on the red links. This is a local operator 
which commutes with all the terms in the Hamiltonian 
except the red shaded plaquette, which sees a flux of 
2 G Z 4 = {0,1, 2,3}. Likewise, we imagine acting with 
ai to the given power on the purple links. This produces 
a charge 2 in the vertex given by the purple diamond. 
Thus we can create even charges and even fluxes locally 
near a Z 2 defect. 


In this subsection we analyze the example of a Z 4 gauge 
theory with G = Z 2 = {1 ,g} acting by p(g) : j -► -j, 
where j G {0,1, 2, 3}. As we discussed earlier, after gaug¬ 
ing G the topological order given by equation [64] is just 
that of the gauge theory of the extension E, where E can 
be either B 8 or (Q) 8 . 

Analyzing extrinsic defects in the ungauged theory is 
the same as including a non-dynamical background G 
gauge field, which just amounts to the following modifi¬ 
cation of the Hamiltonian in equation [64] 


H no„- dy „ am . = _ £ £ Aviv) - £ Bp (67) 
V r]€z%n P 

Note that we are now working with a purely 2d square 
lattice. Let us understand these two terms for a fixed 
background gauge field configuration, namely that of two 
widely separated defects, illustrated in figure [ 6 ] In this 
case, the two terms in the Hamiltonian of equation [67] re¬ 
duce to the ordinary Z 4 gauge theory vertex and plaque¬ 
tte term everywhere except for the plaquettes intersected 
by the dashed blue branch cut line and the vertices di¬ 
rectly above them. For these intersected plaquettes, we 
have a modified plaquette term where ac, is inverted for 
the upper horizontal link in the product over plaquette 
links (recall that is the ‘phase’ operator in the Z 4 
gauge theory). Likewise, just above the dashed blue line 
we have modified vertex terms, where cl (the charge op¬ 
erator in the Z 4 gauge theory) is inverted for the lower 
vertical link (the one bisected by the dashed blue line) in 
the product over vertex links. 

Now, the topological superselection sector of a Z 2 de¬ 
fect can be altered by fusing Z 4 charges and fluxes into 
it. However, as shown in figure [ 6 ] we can create a charge 
of 2 G Z 4 = {0,1, 2, 3} or a flux of 2 G Z 4 using a local 
operator near a Z 2 defect. So the only non-trivial su¬ 
perselection sectors near the Z 2 defect correspond to the 
parity of charges or fluxes bound to the Z 2 defect, and 
there are only 4 such superselection sectors (as opposed 
to the 16 anyon types). This also implies that when we 
fuse a pair of defects, the fusion product can be changed 
by an even charge or even flux with a local operator. 
This implies that a pair of defects can fuse into any even 
charge or even flux, so the fusion rules are non-abelian. 
This is a general feature of defects in permuting theories. 
Indeed, the local operator which creates an even charge in 
the presence of the defect can be thought of as creating 
a pair (fundamental charge, fundamental anti-charge), 
braiding the anti-charge around the defect to turn it into 
a charge, and re-fusing with the leftover charge to make 
a charge of 2 , and similarly for the operator that creates 
an even flux. In the next section we will generalize this 
idea to arbitrary permutation actions. 

Notice that none of the above discussion depended on 
whether we chose E = B 8 or E = Q 8 . We will see in 
the next section that this is a general feature: the su¬ 
perselection sectors of a defect, as well as its fusion rules 
with anyons, depend only on the permutation p, and not 
on w(g,h) or any other data. In order to discriminate 
between E = B 8 and E = Q 8 we have to perform a more 
subtle measurement, involving braiding anyons and de¬ 
fects. We will discuss this in detail in the next section; for 
now just note that this kind of braiding measurement re¬ 
ally comes down to computing commutation relations of 
the string operator that creates defects with other string 
operators, and is not something that can be done in a 
purely fixed background G gauge field configuration. 
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IV. SYMMETRY LOCALIZATION AND 
DEFECT FUSION 

In the previous section we constructed a family of lat¬ 
tice models of SETs where the symmetry G permuted 
the anyons. For each permutation action p, we had mul¬ 
tiple SETs corresponding to different choices of a func¬ 
tion h ) G Hp(G, Z n ). We saw that these SETs were 

distinct by gauging G and examining the full quasipar¬ 
ticle statistics of the gauged theory. However, we also 
saw that the distinctions between these SETs (for fixed 
p but different w(g,h)) were not visible if we examined 
only non-dynamical background G gauge field configura¬ 
tions (that is, with extrinsic defects, but without fluctua¬ 
tions of the G gauge field). In this section, we will study 
the distinctions among such SETs at a formal level, for 
general anyon theories (not necessarily gauge theories) 
and general action of finite on-site unitary symmetry G. 
While we will not fully gauge G, we will consider oper¬ 
ators which create, annihilate, and move defects. This 
algebraic structure turns out to contain the same infor¬ 
mation as the gauged theory, but is somewhat easier to 
deal with, because it is smaller, as it does not contain the 
G charges as topological excitations. 

First, let consider ordinary anyons. The braiding and 
fusion rules of anyons are encoded in a unitary modu¬ 
lar tensor category (UMTC) A, described concretely by 
a collection of quasiparticle labels, fusion spaces, notion 
of antiparticle, and F and R matrices (see appendix E 
of reference f23l for a review). Now, the addition of de¬ 
fects expands the list of labels to the set of all defect- 
anyon composites, and, because defects can also be fused 
and braided, there must be an extension of the algebraic 
structure to this larger set of labels. However, because 
the defects are extrinsic modifications of the Hamiltonian 
rather than excitations, and in particular carry branch 
cuts, this algebraic structure is not a UMTC. Indeed, 
braiding around an extrinsic defect can change anyon 
type and act non-trivially on local operators, which is 
something that cannot happen when we braid around 
an ordinary anyon. The algebraic structure describing 
anyons and defects is a so-called braided G-crossed cat¬ 
egory. See reference \27\ (sec. VI A.-C.) for a detailed 
definition of a braided G-crossed category in the con¬ 
text of defect and anyon braiding. In this paper, we 
will not discuss all of the properties of braided G-crossed 
categories all at once; rather, we will only discuss the 
ones we need, at a physical level, to constrain the pos¬ 
sible defect fusion rules. We build this structure up in 
stages: rather than imposing all of these properties at 
once, we first fix the UMTC of the anyons, then look for 
all consistent solutions for the superselection sectors of a 
single defect and the F an R matrices defining its braid¬ 
ing with anyons, and only then, having fixed this single¬ 
defect structure, study processes involving two defects! 22 ^ 
At each stage, we specify the braided G-crossed proper¬ 
ties that we use. This hierarchical approach parallels the 
mathematical constructions in reference 1201 


In particular, we find that the properties of a single de¬ 
fect, including its topological superselection sectors, and 
braiding and fusion rules with the other anyons, are de¬ 
termined purely by the permutation p (more precisely, by 
the associated action of G by braided auto-equivalences, 
which we define below), and are independent of the frac- 
tionalization class lj. 


A. Non-permuting case 

Let us first give a condensed review of symmetry frac- 
tionalization in the case that G does not permute the 
anyons, i.e. fixes all topological superselection sectors. 
When G does not permute the anyons, we assume the 
existence of an approximate “local” action of G on par¬ 
ticular anyons - see appendix [B] for details on how to con¬ 
struct it. This local action is just an operator U g oc acting 
on the spins within a distance r of the anyon, whose com¬ 
mutation relations with all local operators acting near the 
anyon are approximately the same as those of U g . Here 
the quality of the approximation is controlled by e _r /^, 
where £ is the correlation length; henceforth we take r 
sufficiently large and neglect the exponentially small er¬ 
rors. These properties imply that the U g oc satisfy the 
group relations up to a possible phase ambiguity uj a (g,h) 
(here a is the anyon acted on): 

U l ; c U l ^=u a (g,h)Ull c (68) 

Associativity then gives 

u a (g, h)u a (f, gh) = u a (fg, h)w a {f, g), (69) 

which, together with the uniqueness of the U l ° c up to 
overall phase, imply that we can extract from them a 
unique cohomology class [uj a \ G i^ 2 (G, J7(l)). Further¬ 
more, for every allowed fusion channel a x b c, we must 
have cj a (p, h)c^(p, h) = u c (g,h), since such a fusion is a 
local operation and carries integral symmetry quantum 
numbers (see appendix [B] for a detailed argument). Now 
a technical result of category theory (lemma 3.31 of ref. 
l25l) shows that any J7(l)-valued function of the anyons 
which satisfies this property can be represented as 

h) &u>(g,h),a (70) 

where h ) is an abelian anyon, and Sb, a denotes the 
full braiding phase of the abelian anyon b around a gen¬ 
eral anyon a. From equation [69j we then have that 
[id] G i^ 2 (G, A), where A is the subset of abelian anyons 
viewed as an additive group. Thus we see that the sym¬ 
metry fractionalization data is encoded in a cohomology 
class in [cJ\ G H 2 (G,A). 

Defect fusion rules in the non-permuting case 

Let us now re-interpret this cohomology class [cj] in 
terms of defect fusion rules. As explained in appendix [Bj 
a ^-defect r g is an extrinsic modification of the Hamil¬ 
tonian which inserts a symmetry flux of g in a certain 
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location, with a ^-branch cut emanating from that lo¬ 
cation. We expect a ^-defect to have some topological 
superselection sectors, which can be changed by fusing 
anyons into the core of the defect (in reference f27l these 
are referred to as topologically distinct types of ^-defects, 
and discussed in section V. B.). Now, we will see below 
that in the case of permuting symmetries, there exist 
local processes, involving changing anyon type by braid¬ 
ing around the defect, which seemingly create topological 
charge in the vicinity of the defect; hence, the structure 
of defect superselection sectors in this case depends on 
the permutation being applied. In the non-permuting 
case we are considering now, there are no such processes, 
so we expect that there is no difference between the su¬ 
perselection sectors of r g and those of the trivial defect. 
In other words, the superselection sectors of r g are in one 
to one correspondence with the anyons. 

However, a key point is that there is an ambiguity in 
how to define the Hamiltonian in the core of the ^-defect. 
This means that there is no canonical way to identify de¬ 
fect superselection sectors with anyons: different choices 
of the core Hamiltonian yield different ground state su¬ 
perselection sectors. Now, we can at least demand that 
the core Hamiltonian be chosen in such a way that a pair 
of well separated defects r g and r g - 1 , connected by a 
branch cut, has no ground state degeneracy. This means 
that the ambiguity in the defect ground state superse¬ 
lection sector is at most by an abelian anyon. Having 
no canonical way to resolve this ambiguity, we simply 
choose a ground state superselection sector x g arbitrar¬ 
ily for each ^-defect; this involves making \G\ choices, one 
for each g G G. 

Now consider fusing a ^-defect and an h-defect. This 
is accomplished by a g gauge transformation that merges 
the two branch cuts into a single gh branch cut. The key 
point now is that although this fusion product is a gh de¬ 
fect, because of the ambiguity in the Hamiltonian at the 
core, it may have a different ground state superselection 
sector than the already chosen x g h • Call this difference 
a g ,h', as we just argued, it is an abelian anyon. In fact, 
it is the same as the abelian anyon w(g,h) that was de¬ 
fined above. Indeed, the local action U g oc on some anyon 
b defined above is essentially given by braiding b with a 
^-defect. Then the difference between U g oc Uj^ c and U l g °^ 
is just the difference between braiding b around a pair 
consisting of a g defect and an h defect versus braiding 
b around a single gh defect. According to eq. [68] above, 
this is the braiding phase of b with u(g,h), whereas in 
the present discussion, it is the braiding phase of b with 
a g ,h- Since these two must be equal for all 6, modularity 
implies that a g ^ = h). 

Of course, one has to be a little more precise in defin¬ 
ing the local action of the symmetry: U g oc must be an 
operator that acts on just the anyon 6, and any ^-defects 
involved in the process must be created from the ground 
state and likewise disappear at the end of the process. A 
careful treatment, given in appendix [Bj involves a dual 
picture of holding the anyon b fixed and braiding the 


^-defect around it instead. Nevertheless, the above argu¬ 
ment applies equally well to this definition of U g oc , and 
the conclusion is the same: a g ^ = oj(g, h). 

Thus we can extract the symmetry fractionalization 
data from the defect fusion rules in the non-permuting 
case. 


B. The permuting case: single fixed defect 

The case where G permutes the topological superse¬ 
lection sectors is significantly more complicated. For one 
thing, there is now no canonical ‘trivial’ SET, as op¬ 
posed to the untwisted case, where we can construct a 
trivial SET simply by defining G to act trivially on the 
microscopic degrees of freedom for any Hamiltonian re¬ 
alizing the desired intrinsic topological order. Thus, we 
should not expect to be able to assign group cohomology 
classes to permuting SETs, since there is now no pre¬ 
ferred choice of a trivial SET corresponding to the trivial 
cohomology class. Instead, given a particular choice of 
intrinsic topological order and twisted action of G, our 
approach will be to start with an arbitrary SET realiz¬ 
ing this twisted action, and look for all possible gauge 
inequivalent ways of deforming the defect fusion rules. 
This approach is rooted in the general i dea of c lassify- 
ing SETs via braided G-crossed categories 1 where 
constructing defect fusion rules - formally, a tensor prod¬ 
uct of bimodule categories - is one step in a systematic 
construction of the braided G-crossed category. For now, 
however, let us examine a single ^-defect in the twisted 
case where p(g) non-trivially permutes the anyons. 

To appreciate the complexity of this twisted case, note 
that the ^-defects are now generically non-abelian. For 
example, consider a defect of a Z 2 symmetry which ex¬ 
changes e and m in a model with toric c ode to pological 
order (for an example of such a model see 14 23 ). By nu¬ 
cleating a pair of e’s, braiding one around the defect to 
turn it into an m, and fusing with the remaining e, we 
perform a local process that nucleates the exm = f par¬ 
ticle in the vicinity of the defect. This is reminiscent of an 
Ising defect, and indeed there is a two-fold degeneracy for 
widely separated defects associated with absorbing and 
emitting the / fermion (also, the defects act ually become 
Ising anyons upon gauging the Z 2 symmetr}ff 2 15 27 \ 

To get a handle on this complexity, suppose first that 
we are given a Hamiltonian corresponding to a single 
fixed ^-defect, and suppose further that we know noth¬ 
ing about the location of the branch cut or the action 
of the symmetry on our microscopic degrees of freedom. 
What kind of data can we extract about the symmetry 
enriched phase, using only braiding and fusion operations 
of anyons in the background of this fixed defect configura¬ 
tion? Certainly we can perform experiments that braid 
anyons around the defect and measure the topological 
charge before and after braiding: this is just measuring 
the permutation action of g: 
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FIG. 7: Basic associativity relation for a left module 
category structure of a ^-defect. The ^-defect 
superselection sectors x g ,y g , z g are the simple objects. 
For simplicity we assume in this figure that all fusion 
spaces are one dimensional. The sum on the right hand 
side is over all anyons c which satisfy the fusion rules, 
i.e. correspond to non-trivial splitting spaces V^ ,b and 

T y c i x g 

Vz 9 • 


a —y g • a. (71) 

Furthermore, we can also measure the action of g on 
anyon fusion spaces: 

g : V?’ b -> Vg 9 c’ 9 ' b (72) 

X eV c a ’ b ^g-X e V g 9 c a ’ 9 ' b , (73) 

A physical process that can measure this action on anyon 
fusion spaces is as follows. Given an operator X which 
splits c into the pair a, 6, i.e. X G V^’ 6 , define g • X be 
the operator that acts as follows: first braid g • a and 
g • b around the defect to turn them into a, b respectively, 
apply X, and then braid the fusion product c in the op¬ 
posite direction to get g-c. Note that because braiding an 
anyon a around the defect is only defined up to a phase 
Ota-, g ’ X is only well defined up to a phase a^a^a^ 1 - 
A permutation action on anyons together with a com¬ 
patible action by unitary linear transformations on anyon 
fusion spaces is called a braided tensor autoequivalence of 
A (in reference|27, this also goes under the name of topo¬ 
logical symmetry, discussed in section III. A.). There is 
also a notion of two braided autoequivalences being the 
same: this occurs when the two autoequivalences have 
the same permutation action, and their actions on the 
fusion space V^ ,b differ by a^a^a^ 1 , where the {a a } is 
some fixed set of 17(1) phases. Hence, by performing 
braiding experiments around a fixed g defect configura¬ 
tion we uniquely recover precisely a braided tensor au- 
toequivalence corresponding to g , up to this notion of 
braided tensor autoequivalences being the sameP^ 

Now, the goal of the classification program of SETs is 
to construct a larger algebraic structure - the braided G- 
crossed category - which describes the fusion and braiding 
rules of both anyons and defects. How much of this larger 
structure can one recover from knowing the braided ten¬ 
sor autoequivalence corresponding to each g G G? The 
answer is given by theorem 5.2 of reference [20} the 
braided tensor autoequivalence uniquely determines an 


‘invertible bimodule category’ corresponding to the g- 
defect. Roughly, the ‘invertible bimodule category’ cor¬ 
responding to the ^-defect is the subset of the braiding 
and fusion data of the theory that involves a single g- 
defect and an arbitrary number of anyons, modulo gauge 
equivalences. The intuition is that, while it may be pos¬ 
sible to find a gauge equivalence between such a subset 
of braiding and fusion data for two different theories, it 
might be impossible to extend this gauge equivalence to 
a gauge equivalence of the full theories, involving braid¬ 
ing and fusion data of an arbitrary number of defects and 
anyons. So the invertible bimodule categories associated 
to the various defects constitute coarse data that only 
partially classifies different SETs. In particular, in the 
context of the lattice models from the previous section, 
this data is sensitive only to the permutation p, and not 
to the co-cycle uj. 

Let us now define invertible bimodule categories in a 
little more detail. 

Invertible bimodule categories 

First of all, let us put all of our defects and anyons on 
a line. Then we can hold the defect at a fixed position, 
and imagine fusing anyons into and out of it. We can 
do this either from the left or the right; the correspond¬ 
ing structures, consisting of the superselection sectors of 
the p-defect, together with the F-move isomorphisms for 
these fusion processes, are called left and right module 
categories respectively. More precisely, a (say) left mod¬ 
ule category involves this fusion and F-move data for 
anyons fusing in from the left, modulo gauge equivalence 
corresponding to basis redefinitions within the anyon- 
defect fusion spaces, and satisfying pentagon equation 
constraints involving 3 anyons and one defect. The data 
defining a left module category is shown in figure [7] 

Now, physically we can extract the data of both a left 
and right module category structure for a single fixed de¬ 
fect. However, in a physical system these data are not 
independent: there are additional data and compatibil¬ 
ity constraints between fusion from the left and right, 
beyond those encoded in the left and right module struc¬ 
tures alone. This is the associativity of a process where 
one anyon fuses in from the left and another fuses in from 
the right, constrained by the appropriate pentagon equa¬ 
tions. A left and right module category together with 
such additional data and constraints is called a bimodule 
category. Thus a single defect gives rise to a bimodule 
category over the anyons. 

Additionally though, all of the bimodule categories 
that we encounter have further properties due to the exis¬ 
tence of a braiding structure. For now this braiding struc¬ 
ture involves leaving the p-defect fixed, and just braiding 
anyons around it. Let us imagine that the branch cut 
emanating from the defect goes towards the back of the 
page (away from the line where all the quasiparticles sit). 
Then braiding an anyon a in front of the defect, so as to 
avoid its branch cut, gives an isomorphism of anyon de¬ 
fect fusion spaces: 
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TDCL,Xg . ya,x g _v yx g ,a 

y 9 • v y g ^ v y g 


(74) 


Here x g and y g are defect superselection sectors and a is 
an anyon. We can identify the two fusion spaces V g ^ 9 


and Vyg 9 ' using the isomorphism in eq. 74, as illustrated 
in figure [ 8 | Then the compatibility of oraiding and fu¬ 
sion, formally expressed in terms of the hexagon equa¬ 
tions, uniquely determines all the F-move isomorphisms 
involving two anyons and one defect in terms of those 
where the anyons only fuse in from the left, illustrated 
in figure [7] Mathematically one says that in this case 
the bimodule category structure is induced from a left 
module category via the braiding in A. Thus a single 
defect gives rise to a bimodule category over the anyons, 
where the bimodule structure is induced form a left mod¬ 
ule structure by using the anyon braiding. 

A second property of these bimodule categories comes 
from the fact that anyons can braid behind the defect 
as well. Various hexagon equation constraints^ turn out 
to determine the half-braiding of an anyon behind a de¬ 
fect up to a U{1 ) phase which depends on the anyon be¬ 
ing braided, but not on the superselection sector of the 
defect®! These are just the phases a a mentioned ear¬ 
lier. Mathematically, these extra constraints imply that 
the bimodule category is ‘invertible’ - see reference [ 20 ] 
for a precise definition of invertibility, which we will not 
state here. Physically, however, invertibility turns out to 
mean that for any defect, there exists an anti-defect that 
can annihilate with it, under the defect fusion product 
defined in the next section. 

Thus, in summary, just from braiding experiments in¬ 
volving a single static ^-defect one can determine the 
action of g by braided tensor autoequivalences, which 
is equivalent to all fusion and braiding rules involving 
a single ^-defect, with the exception of braiding behind 
the defect, which is determined only up to a phase that 
depends on the anyon being braided (but not the su¬ 
perselection sectors of the defect). The uniqueness here 
is up to a gauge freedom corresponding to a change of 
orthonormal bases in the anyon-defect and defect-anyon 
fusion spaces. Mathematically, this collection of data is 
an invertible bimodule category with bimodule structure 
induced from braiding in the anyons A. 


Example: Z4 gauge theory with Z 2 symmetry 

Let us give a non-trivial example of the above struc¬ 
ture. Let G = Z 2 , and consider Z4 gauge theory, whose 
anyon content is Z4 x Z4. This is the theory we discussed 
in the previous section. We will think of the first Z4 as 
charge and the second Z4 as flux, and denote anyons by 
0 <i,j <4. We will take the non-trivial generator 
g G Z 2 to act on the anyons by (i, j ) (4—i, 4 —j). What 

is the invertible bimodule category structure associated 
to a ^-defect? 

First of all, by the argument above, we can nucle¬ 
ate anyons of the form (2i,2j) in the presence of the 
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FIG. 8 : In the algebraic fusion and braiding theory of 
anyons and defects, one considers two separate splitting 
spaces Vyg Xa and V^ 9 ’ a , corresponding to the anyon a 
splitting off from the defect either to the left or to the 
right. Because of the existence of braiding, these spaces 
are actually isomorphic: Vy g Xg = Vy g 9,a . In our 
graphical calculus, we fix this isomorphism to be the 
one given by braiding the anyon a in front of the defect, 
so as to avoid the defect branch cut, which extends into 
the page. Having chosen a basis for the same anyon 
splitting from the left, this fixes the basis for the fusion 
space of the anyon splitting from the right, one in which 
the half-braid R-move is equal to the identity: 


^-defect. Thus the ^-defect has only 4 superselection sec¬ 
tors, namely (a,/3) g , 0 < a,/3 < 1. The fusion rules with 
the anyons are simply given by addition mod 2 

(i, j) x (a, /3) g = {a + i mod 2 , /3 + j mod 2 ) g . (75) 

These are the same fusion rules that we found directly in 
the previous section. 

The F symbols involving one defect and two anyons 
are non-trivial. There is a choice of basis in which they 
can all be set equal to ±1, and they can be determined 
explicitly using a general equivalence between such in¬ 
vertible bimodule categories and Lagrangian subgroups 
of i/0 H* for any abelian H gauge theory - see the dis¬ 
cussion in section 10.3 of reference uni More explicitly, 
in a particular gauge they can be taken as follows. Given 
even anyons a\ = ( 2 ii, 2 j 1 ), a 2 = ( 2 i 2 , 2 j 2 ), let 

pa!,a 2 ,x _ (_iyij2-t2ji {76) 

Anyons not of this form can always be uniquely written 
as a\ = ( 2 ii+/ii, 2 jq-fz/i),a 2 = ( 2 i 2 +Mi> 2 j 2 + ^i), where 
Mi, 2 , Ml ,2 = 0,1. For these anyons we define the F-symbol 
the same way, Fy 1,CL2,x = (— \y^ 2 -i 2 ji m F-matrices for 
fusion of anyons from the right are then uniquely deter¬ 
mined using the rule in figure [ 8 ] and ordinary anyon fusion 
rules. The R-move corresponding to braiding an anyon a 
behind a defect is also determined uniquely up to a phase 
that depends on a but not on the superselection sector 
of the defect, as is the action on the anyon fusion spaces, 
which turns out to be trivial. 
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Let us emphasize that the F-matrices we have written 
down here, involving one defect and two anyons, will be 
valid for both of the SETs with this permutation action 
on the anyons. The data which distinguishes between 
the two SETs, namely the cohomology classes c o±(g,h) 
discussed in the previous section, will show up only in 
quantities involving two defects. Below we will discuss 
this quantity uj in our more general context. 


C. Defect fusion in the twisted case 


Above we saw that the gauge invariant data contained 
in a single ^-defect is entirely determined by the action of 
g by braided autoequivalences on the anyons. In partic¬ 
ular, we cannot distinguish among SETs with the same 
braided autoequivalence action of G by considering only 
single defects. In this section, we will see that such SETs 
can be distinguished by considering pairs of defects and 
their splitting and fusion rules. 

First let us fix a pair g,h E G, and consider a gh defect 
splitting into a g defect and an ft defect. One can again 
look at associativity relations for processes which start 
with the gh defect and end in a g defect, an h defect, 
and an additional anyon a. There are three types of 
processes, depending on whether a appears to the left, 
to the right, or in between the defects in our chosen line 
ordering; they are illustrated in figure [9| 



FIG. 9: Data defining a defect fusion product of a g 
defect and an h defect into a gh defect. Dotted lines 
represent anyons, solid lines represent defects. 


The coherence conditions on the F-matrices illustrated 
in figure [9] are highly constraining. Indeed, these F- 
matrices must satisfy compatibility conditions with each 


other, with the invertible module category structure, and 
with the braiding and fusion structure of the original 
anyons. These all take the form of pentagon equations 
and variants of the hexagon equations involving 0, 1, or 
2 defects. Despite the complicated nature of these equa¬ 
tions, reference HOj again proves a theorem that constrains 
the possible solutions. Specifically, proposition 7.3(h) of 
l20lshows that, for a fixed pair (g,ft), and having fixed all 
of the one defect data (namely, the map from G into the 
braided autoequivalences of A), the set of gauge equiva¬ 
lence classes of solutions for the fusion rules of g and h 
defects into a gh defect together with the associated F- 
matrices is in one to one correspondence with the set of 
abelian anyons A. Indeed, given one solution, reference 
20 shows how to ‘deform’ it by any abelian anyon a g ^. 
Intuitively, one can think of this deformation as splitting 
off a Q: h from the gh defect before it splits into g and h - 
see figure [lOj We stress that this intuitive picture is just 
graphical shorthand for a precise formula, written down 
in appendix [Cj that expresses the new, deformed fusion 
rules and F-matrices in terms of the original ones. 
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FIG. 10: Given a valid defect fusion product, one can 
construct a deformed defect fusion product, shown here 
with a tilde, as shown in the above figure. The equation 
in the above figure amounts to a construction of new 
F-symbols involving the g and h defects and arbitrary 
anyons. For details on how these new F-symbols are 
defined, see appendix |C| 


Now, suppose we have a complete solution for all of 
the data in our theory, including F-matrices with any 
number of anyons or defects. We can try to find a new 
solution by leaving alone all data involving at most one 
defect (in particular the fusion and braiding data of the 
anyons and the invertible bimodule category structure), 
deforming, for each pair (g,ft), the fusion and F-matrix 
data involving two defects g and ft by some anyon a^, 
as in figure [lO] and attempting to solve for the rest of 
the structure (including associativity conditions with 3 
defects, and the defect braiding matrices). That is, we 
perform the deformation discussed above for every pair 
of defects g and ft. Now, there is no guarantee that if we 
do this then we will be able to solve all of the remaining 
coherence conditions and obtain a consistent theory. A 
basic requirement on a g ^ is 

Pf( a gih)Uf,gh = a fg,h&f,g (77) 

where pf represents the permutation action of / on the 
anyon labels. In the non-permuting case this condition is 
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simply the requirement that we obtain the same topo¬ 
logical superselection sector regardless of the order in 
which we fuse three defects /, g,h. In the twisted case, 
equation 77 is not as obvious, but rather follows from 
imposing all pentagon equations involving three defects 
and one anyon: a solution for the three defect F-matrix 
can be found only if eq. [77] is satisfied. In the language 
of reference [ 20 ] these pentagon equations translate to the 
statement that the tensor product of bimodule categories 
referred to previously is associative, and the condition 
follows from their Theorem 8.5. We note that even when 
this is so, this associativity isomorphism might not satisfy 
the 4-defect pentagon equation and fail to give a consis¬ 
tent theory 32 - however, we will not be concerned with 
this issue here. 


Because deformations of the form a g ^ = Pg(bh)b~^b g , 
where the b g are abelian, can be undone by a re-labeling 
of defect superselection sectors, eq. [77] shows that the po¬ 
tentially non-trivial deformations correspond to nonzero 
twisted cohomology classes in [cj] G i7 2 wisted (G, ^l), rep¬ 
resented by the functions w(g, h ) = a g ,h- In section IV of 
reference there is defined a collection of [/(l)-valued 
2 -cocycles r a (g, h), where a ranges over all of the anyons 
in the theory (not just abelian ones). This r a (g,h) is 
related to our u(g, h ) by: 


ZaG?? &a,co(g,h ) (^8) 


where S a ^ g ^ is the full U( 1) braiding phase of the two 
anyons. 

Of course, given such a nonzero cohomology class, it is 
not a priori clear how many consistent theories it leads to, 
if any. Nevertheless, assuming there does exist at least 
one consistent deformed solution for a given [cj] , one can 
ask what gauge invariant quantity distinguishes this de¬ 
formed solution from the un-deformed one. In principle, 
one can simply compute all of the relevant (two defect, 
one anyon) F-matrices for both solutions and check if 
there is a gauge equivalence between them. However, in 
practice this is difficult, both because of the large number 
of F-matrices involved, and also because extracting uni¬ 
versal data associated to F-matrices is significantly more 
subtle than extracting braiding data for a given physical 
system. 

Fortunately, it turns out that once we specify the defor¬ 
mation [cj], various coherence conditions - namely pen¬ 
tagon and analogues of the hexagon equations - force 
much of the remaining fusion and braiding structure. In¬ 
deed, according to the general classification of references 
l26l and the SET is then actually determined almost 
uniquely: the only ambiguity is the stacking of an ex¬ 
tra G-SPT (symmetry protected topological phase) on 
top of the system. In particular, the braiding of anyons 
around defects is completely determined. As we will see 
in the example below, the gauge invariants that can be 
constructed from this braiding data contain a much more 
convenient signature of the deformation. 


D. Example: Z4 gauge theory with Z2 symmetry 

Let us determine the defect fusion rules in the example 
discussed above. The only non-trivial case is a ^-defect 
fusing with itself, and in this case, since we know that 
braiding around two such ^-defect induces the trivial per¬ 
mutation on anyon labels and the trivial map on anyon 
fusion spaces, from theorem 5.2 of reference ’20 it follows 
that the product of two ^-defects is a trivial defect. Re¬ 
call that we are representing the anyons by pairs (i, j), 
where i denotes charge Z 4 charge and j denotes Z 4 flux, 
that the superselection sectors of the non-trivial ^-defect 
are labeled by pairs (a,/3) p , with a,/3 defined modulo 2 , 
and that the fusion rules are: 


(a,/3) g x (a',P') g = j), (79) 

where the sum on the right is over 0 < i,j < 4 with 
a + a' = i mod 2 and (3 + /3' = j mod 2 . Notice that these 
fusion rules are non-abelian. 

Now, according to the above discussion, the potential 
deformations of this defect product are parametrized by 
^twisted0^2 , ^4 x ^ 4 ) = ^2 X ^2- Let us now focus on 
the deformations where the anyons a g ^ are fluxes of the 
Z 4 gauge theory, which restricts us to a Z 2 subgroup of 
Z 2 xZ 2 ; the non-trivial co-cycle is then a g , g = (0, 2 ). We 
do this because this Z 2 is the one that can be accessed in 
the particular class of exactly solved models given in the 
previous section. Indeed, these models, like those of ref¬ 
erence [l]on which they are based, were constructed to al¬ 
low for fractional symmetry charge assignments to gauge 
theory charges, since it is the gauge charges that are sen¬ 
sitive to the Berry flux in the vertical plaquette terms. 
The other Z 2 would correspond to a phase where instead 
the gauge fluxes carry fractional symmetry charge, and 
would require a different lattice construction 29 . 

Note that deforming by this co-cycle does not change 
the fusion rules in equation [79] Thus it is now much 
more difficult to distinguish the deformed and unde¬ 
formed products of bimodule categories. Indeed, as op¬ 
posed to the untwisted case, where the distinction can 
be ascertained by just examining the fusion rules for de¬ 
fect superselection sectors, we now have to examine sub¬ 
tle gauge invariant F-matrix data for processes involving 
two defects and one anyon. Fortunately, as discussed 
above, this distinction appears also in the braiding data 
of anyons with the defect in the full solutions (‘braided 
G-crossed categories’) corresponding to the deformed and 
undeformed data. As shown in appendix [Cj imposing 
hexagon coherence conditions shows that the deformed 
and undeformed cases differ in the sign of the gauge in¬ 
variant braiding quantity 


Rxg,(3,0)R(3,0),XgRx gf (l,0)R(l,0),Xg, (^ 0 ) 

where x g is an arbitrary superselection sector of the g- 
defect, and R a ^ is the R-matrix corresponding to the 
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exchange of a and b. Recall that (i,j) represents an 
anyon which is a composite of a charge i and flux j 
i,j = 0,... 5 3. Thus equation |80| simply gives the braid¬ 
ing phase for the process of braiding a fundamental Z 4 
charge (1,0) twice around the ^-defect. Notice that 
braiding it only once turns it into the opposite charge 
(3,0) and hence does not lead to a gauge invariant phase. 


Now recall that in the previous section we had an ex¬ 
actly solved lattice construction for precisely this topo¬ 
logical order and action of symmetry, and found two dis¬ 
tinct SETs. We claim that these are lattice realizations of 
the two cases - deformed and un-deformed - that we are 
discussing right now. In order to argue this, we will just 
show that the two distinct lattice SETs we constructed 
in the previous section have a different sign of the gauge 
invariant quantity defined in eq. [80] Recall that upon 
gauging G, the two SETs we constructed in the previous 
section resulted in B 4 and Qg topological orders, respec¬ 
tively. 


First, let us establish some more notation. The second 
group cohomology H 2 {^ 2 ^ 4 ) = Z2, where Z2 = {1 , g} 
acts as p : x — >• x~ x on Z4 = {1, z, —1, — i}. We found 
two representative co-cycles, 0 ;+ and c c was equal 
to 1 for all values of its arguments, and uj-(g,g) = — 1 . 
According to Appendix [Dj uj+ corresponds to the di¬ 
hedral group B4, and Co>_ to the quaternion group Qg. 
The first of these just consists of the symmetries of the 
square, with a Z4 rotation subgroup which we refer to as 
Z4 = — 1, — i}, while the second consists of the el¬ 

ements {± 1 , =bi, ± 7 , =b/c} which follow the multiplication 
rules of the quaternion group. Here we again take the Z 4 
subgroup to be ±1, ±i. The key difference between these 
two extensions, captured by the respective co-cycles, is 
that in B4, any lift of g squares to +1, whereas in Qg, 
any such lift squares to — 1 . 


We now perform the double braid computation in eq. 
|80| in the fully gauged theories. Here, the charges (1,0) 
and (3,0) (obtained from it by acting with g) pair up 
into the unique two dimensional irreducible representa¬ 
tion of either B 4 or Qg depending on whether we are in 
the l or case, respectively. Now, the gauge theories 
(or ‘quantum doubles’) of these two groups can be con¬ 
structed using a general prescription outlined e.g. in sec¬ 
tion 3.2 of reference l33l The key formula for computing 
braiding is the general formula for the universal R ma¬ 
trix given there. Here we are interested in moving a pure 
charge - our two dimensional irreducible representation 
of B 4 or Qg - around a g-flux twice. Using the formula 
for the R matrix given in section 3.2 of reference [33j the 
action of such a double braid operation is simply given by 
g 2 acting within the two dimensional irreducible repre¬ 
sentation of E that contains the fundamental Z4 charge. 
This is just the identity in B 4 , and —1 in Qg, which acts 
as —1 in the two dimensional irreducible representation, 
so that the value of the double braid in eq. 80 differs by 
a sign for the two theories, as desired. 


V. SUMMARY AND NEW DIRECTIONS 

In this paper we presented a framework for studying 2d 
bosonic symmetry enriched topological phases with a fi¬ 
nite, unitary symmetry group G. We demonstrated that 
even when G acts non-trivially on the topological supers¬ 
election sectors, we can still extract fractionalization data 
and distinguish different symmetry enriched phases using 
the fusion rules of extrinsic twist defects. Moreover, some 
non-trivial features such fractionalization data are visible 
only at the level of associativity rules - i.e. F-move iso¬ 
morphisms. We illustrated these results in a class of ex¬ 
actly solved lattice models of Z n gauge theory, for which 
the fractionalization data, in the form of a twisted group 
cohomology class [cj] £ H 2 (G, Z n ), is one piece of input 
defining the Hamiltonian. By gauging G we found that 
the intrinsic topological order of the resulting model is 
that of an E gauge theory, where E is the extension of G 
by H determined by the permutation action p and [uj]. 
In particular, whenever these gauge theories are distinct, 
so are the underlying SETs. 

In the future we would like to extend our results to the 
case of time reversal symmetry. Its anti-unitary nature 
makes it difficult to define the notion of a twist defect 
and of gauging, but there are indications 34 that some 
notion of symmetry fractionalization should still hold in 
this case. In particular, we believe that the exactly- 
solved models of reference |T| can be generalized to this 
case. Also, we would like to study 3 dimensional gen¬ 
eralizations of these models. Indeed, topological order 
and symmetry enrichment are not as well understood in 
3 dimensions as they are in 2 dimensions, and it would 
be nice to have a class of exactly solved 3d models to 
benchmark potential classification schemes. Another set 
of systems we would like to understand are the Walker- 
Wang models of 3d SPTs - perhaps the ideas presented 
in this paper can be used to give a more geometric con¬ 
struction of the symmetry decorated semion model of 
reference [32, Finally, it would be good to see whether 
these ideas generalize to fermionic systems, and in par¬ 
ticular, to what extent it is useful to think of gauging the 
global fermion parity symmetry. 
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Appendix A: Derivation of E gauge theory 
Hamiltonian 


In this sub-section we will determine the low energy 


57 


physics of the gauged Hamiltonian given in eq. 
Again, this sub-section is rather technical, but the fi¬ 
nal result is quite simple: the low energy physics of this 
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gauged Hamiltonian is equivalent to that of the non- 
abelian lattice gauge theory based on the extension E 
of G by Z n determined by p, given in equation |A22| 

Our analysis here is similar to that of reference [T| in 
the non-anyon-permuting case. The general strategy is 
to simplify the problem by imposing certain terms in 
the gauged Hamiltonian in eq. [57] as constraints. This 
gets rid of states in the Hilbert space which violate these 
constraints, leaving us with a smaller Hilbert space that 
is hopefully easier to understand. However, we have to 
be careful and ensure that by doing this we have not 
removed any topologically non-trivial excitations from 
the spectrum, because if we have, we no longer have a 
model with the same topological order. A second issue 
is that the resulting Hilbert space, although smaller and 
in principle easier to understand than the original, is a 
constrained Hilbert space. That is, it is not of the form 
of a tensor product of site Hilbert spaces. Because we 
would like to work with unconstrained Hilbert spaces, 
we resolve this issue by rewriting the constrained Hilbert 
space as an unconstrained Hilbert space together with 
a reparametrization redundancy. Again, we show that 
violating this reparametrization redundancy does not in¬ 
troduce additional topological superselection sectors, so 
that that the reparametrization redundancy condition 
can be imposed energetically. Thus, at the end we have 
another generalized spin model with the same topological 
order and low energy physics; in effect we have carried 
out a duality transformation. This new generalized spin 
model, after some trivial manipulations, turns out to be 
the gauge theory of the extension E. 


1. Establishing the Hilbert space constraints 

a. Supervertex constraints 


We first impose as constraints the plaquette terms C p 
for plaquettes p contained entirely within a single super¬ 
vertex V: that is, we restrict to the states |T) for which 
C p |T) = |T). For each plaquette triple (/, g, h), this con¬ 
straint reads 

w(f,g)w(fg,h)ri V (f,gh)~ 1 = u)(g,h) p{f) (Al) 

Using eq. [32j we see that one solution to this equa tion 
is fj(g,h) = uo(g,h). Now, any other solution to eq. |A1 
must be related to p(g, h ) by a gauge transformation: 


Vv(d, h) = 0 V ( gh)rj(g , h)0 v (g) 


(A2) 


We now define an unconstrained Hilbert space by trad¬ 
ing in the Z n link variables rjv(g,h), i.e. those with 
links entirely within supervertex V, for Z n valued ‘clock’ 
variables 6y{g )• We introduce the usual operator alge¬ 
bra associated to each Z n clock: Oy(g) is the opera¬ 
tor that has 0y(g) as an eigenvalue, and <&v(g) acts as 
I Ov(g')) -A |e “ 9 ' Oy(g')). One advantage of r epresent- 

ing gy(g , h) in terms of the 0y(g) as in eq. A2 is that it 


simplifies the vertex term Ay(g). Indeed, the contribu¬ 
tions to Ay(g) coming from links within the supervertex 
V amount to just ^y 1 ^), so the vertex term becomes: 

A v {g) -> $y(sO tl e L {g) s ^^ L) (A3) 

L~y 


Note that { 0y{g )} form a redundant parametrization of 
the original link degrees of freedom, in that shifting all 
0 y(g) simultaneously by the same ^-independent amount 
does not change any rjy(g , h) defined in eq. A2 Thus the 
Hilbert space of the gauged Hamiltonian in eq. [57] with 
the supervertex constraints C p |T) = |T),pGU imposed, 
is recovered by projecting onto the subspace where 


n*?GOI¥>H¥>,v m ez„ (A4) 

geo 


This is the reparametrization redundancy condition we 
mentioned earlier. 

What about the Hamiltonian? We can certainly write 


down a Hamiltonian which reduces to that in eq. 57 


when acting on the subspace satisfying the condition in 
eq. |A4[ and which also includes projectors onto this sub¬ 
space, such that these projectors commute with the rest 
of the Hamiltonian. However, before we do this we will 
impose the second set of constraints in the next subsec¬ 
tion. For now, however, we just note that such a Hamil¬ 
tonian will necessarily have the same topological order 
as the original one in eq. 50 To see this, we argue as 
follows. First of all, one may worry that there are non¬ 
trivial topological excitations in the Hamiltonian of eq. 


50 which are removed from the spectrum once the con¬ 


straint C p = 1 for p G V is imposed. However, it is 
fairly easy to see that this is not the case. Indeed, any 
violation of C p associated to V can be removed with a 
local operator supported only on V. This operator may 
create other excitations which are violations of C p > for 
p' containing some horizontal links, but crucially, it will 
not create violations of any C p " for p" contained entirely 
in some supervertex. Now, the violations of C p > for p' 
containing horizontal links can again be removed by act¬ 
ing on the appropriate horizontal link variables. These 
will not create excitations at any other C p . Thus, any 
local excitation can be pushed off to violating only the 
B p (g) and remaining terms in the Hamiltonian using a 
local process, so that its topological superselection sec¬ 
tor remains in the spectrum even after the constraint is 
imposed. A similar argument shows that imposing the 
reparametrization redundancy condition in eq. |A4| does 
not introduce new topological super select ion sectors. 


b. Superlink constraints 


Let us now also impose the constraints C p |T) = |T) 
for p involving a superlink L = (VV'). The definition of 
Cp in eq. 53 shows that C v just measures the Z n flux in a 
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certain rectangular plaquette involving the superlink L. 
A slight rearrangement of eq. [53] shows that setting this 
Z n flux to 1 is equivalent to: 

VL{gh)ri v {gl l g,h) p(9L) = w(g,h)r] L (g) 


(see eq. |A2| ) we obtain: 

V , L(g h ,9L)v{gZ 1 g,h) p{9L) = f}(g,h)rj' L (g,g L ) 
where rf L (g,g L ) = 0 V '{g)riL{g)0v{gZ l g)~ p(9L) (A7) 


one to fix is rj f L (gL, Ql), simplifying eq. A6 


to 


^(SzAflx) =u(g L ,h)r]' L (gL,gL) 


From this point forward, we will refer to r}' L (gL,9L ) as 
r ]' L , i.e. drop the two arguments. Imposing the super¬ 
link plaquette constraint thus leads to a Hilbert space 
which can be parametrized by introducing a single Z n 
clock variable g' L for each superlink L, in place of the 
link variables rjL(g)- Note that there is no additional 
gauge redundancy introduced by this parametrization. 
We let a' L and e' L be the evaluation and shift operators 
corresponding to the g' L variable. 


9l 



9 h 


,9 


FIG. 11: Links acted upon by the superlink constraint 


2. Hamiltonian after imposing constraints on the 
gauged theory 

We have thus constructed a Hilbert space parametrized 
by {Ov(g)}v, Wli9l}l- Our goal here is to write out the 
Hamiltonian in this Hilbert space. We begin with the 
plaquette term Bp(g ) in eq. 


52 


assuming, without loss 
of generality, that all links are oriented positively relative 
to the plaquette p, which we take to be j-sided. Order 
the j links 1 ,... ,j starting from an arbitrary point, and 
for convenience let gi = denote the G gauge field 


variable on the ith link, 1 < i < j. In appendix 1 we 
show that the following term reproduces Bp(g) in the 
new Hilbert space: 


(A5) 



terms of 

Tr G,z n B p (g) -> Tr z n ^ P B' p Tr gB p 

(A9) 

s section 

_ ij \ 



^p = 9ri 9i J 

(A10) 


i=l \r= 1 / 


(A6) 
b (A7) 


(AH) 

g L ) for a 

3 


A useful 

Bp = JJ Oti 

(A12) 

(AS) 

where Tr^ refers to tracing over the regular representa- 

tion of G. 



Now, for each superlink L, the pair ( 77 ^, #l) can nat¬ 
urally be thought of as an element of the extension E 
of G by Z n determined by cj(p, h), since the extension is 
just a twisted product of Z n and G (see appendix [P] for 
some background on group extensions). The twisting in 
this product is manifested only in the non-trivial nature 
of the group multiplication law in E\ 


(v' 1,91 ) • 0?2>52) = ( vW 2 Pi9l) u(gi,g 2 ),gig 2 ) (A13) 

This way of thinking about the degrees of freedom on 
superlinks is especially convenient because the plaque¬ 
tte term in eq. |A9| precisely just enforces that the pairs 
around a plaquette multiply to the identity in 
the extension E. Indeed, repeated applications of the 
above group multiplication law in E give 


ri(^) = if^ 1 ^ n 


n gr,gi j ,in) 

Vr=l / i / 

(A14) 


which is precisely the quantity that the term in eq. |A9 
energetically prefers to be trivial. Thus we just have the 
plaquette term of a non-abelian E gauge theory. 

As for the vertex term, a quick calculation shows the 
action of Ay(g) within the constrained Hilbert space is 
the same as that of the operator &' v (g) defined by: 


&v(g) I 0 I °v{g'))® 0 Wl,9l) | = 


\e ™ 5 °'<i'0 v (g')) 


I Vl,9l) (A15) 


g'eG 


L~V 


This is a great simplification, because &' v (g) does not act 
at all on the g' L degrees of freedom. On the other hand, 
the constraint |A4[ which now also has to be imposed 
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energetically, is less trivial. We have: 


(8) 1 9v{g')) ® (^) I Vl>9l) 

geo \g'€G 


L~V 


0|^(5)e^)®0|e 
. 01 0 L ( 9 )e™)®®W L e 


Vl’9l),L 


V 

,9l),L<- V 

(A16) 


This is equal to the action of the operator 


A'v ®v(g) 

g£G 


with 


where 


Ay = (e f L ) 

L~V 


r{L) 


c'lWl>9l) = I e 2 *n' L ,g L ) 

(' s'l) R \Vl>9l ) = ,9 l) 


(A17) 

(A18) 

(A19) 

(A 20 ) 


Now, since Yl g &y(g) already appears in the Hamilto¬ 
nian, and since all the terms in the Hamiltonian will 
commute, we can replace it with 1 in the expression 
AyYlgeG^vid)’ so ultimately the constraint 


A4 


posed by including the terms— V m Ay 1 in the Hamilto¬ 


nian. Now, according to eq. A18, A'™ is just the gauge 
transformation by m G Z n C E at supervertex V, when 
we interpret the degrees of freedom as those of an E gauge 
theory. Thus so far we have recovered the plaquette term 
of the E gauge theory, and a term that imposes gauge 
invariance with respect to the Z n subgroup of E. 

Finally, we must write down terms corresponding to 
Ay(g), which make the G gauge field fluctuate. There 
is some ambiguity in the choice of these terms because 
the action on the portion of the Hilbert space violating 
the constraint | A4| is not uniquely determined. However, 
one valid choice, which commutes with this constraint 
and the rest of the terms in the Hamiltonian, is to sim¬ 
ply define Ay{g) to act as the gauge transformation by 
(1,#) G E ; the arbitrariness is due to the fact that we 
could have equally chosen (e 27Tirn / n , g) G E. 

So, our Hamiltonian, when written in terms of opera¬ 
tors that act on {Oy (g), 77 ^, gq,}, becomes 


= - E E (<?) - E Tr Tr gBp 

V m£Z n P 

-EE A 'v-EE Ay(g) (A 21 ) 

v m y g£G 

This is almost the Hamiltonian of an E gauge theory. 
Indeed, the first term simply enforces that the ground 
state is symmetrized over the E s, and the 6 degrees of 
freedom do not appear anywhere else in the Hamiltonian. 
Thus we can remove them from the Hamiltonian without 


changing the low energy physics or the topological order. 
As far as the rest of the Hamiltonian, we can make it 
exactly equal to that of an E gauge theory by adding in 
terms of the form A , y l Ay{g ); these do not introduce any 
new conditions: 


= Avia) 

v gee 

m£Z n 

- Tr Z „ VpB'p Tr G Bp (A22) 

P 

This is simply the Hamiltonian of the discrete lattice 
gauge theory based on the group E, the extension of G 
by Z n defined by our group cohomology class u(g, h). In 
particular, A' v and Ay(g) satisfy 

A v {gMv{92) = A'f 9l ’ 92) Av(gi92 ) (A23) 

where e 2 % lq ( 91 ’ 92 ') = u;(< 71 , 02 ) 
A v (g)Ay = A'f 9)m A v (g) (A24) 

The Hamiltonian H# is equivalent to that in eq. [64] 
While the arguments we have given were written down 
only for Z n gauge theory, it is quite easy - though nota- 
tionally cumbersome - to generalize to the gauge theory 
of any finite abelian group. 


Appendix B: Symmetry fractionalization, local 
action of G , and defect fusion rules 

1 . Local action of G on anyons 

Throughout this subsection, we assume that G fixes 
the topological superselection sectors. Suppose we have 
a state with anyon a at position x, and its anti-particle 
somewhere far away at position y. In order to extract 
the fractional quantum numbers, we would like to define 
a local action of g G G in a neighborhood of x - not in¬ 
cluding y - so that we are sensitive to the Berry phases 
associated with acting on a, without also picking up the 
opposite phases from a. More precisely, for each g G G 
we would like to have an operator U g that acts only on 
spins near x, that commutes with the Hamiltonian, and 
that has the same commutation relations with operators 
localized near x as the global action U g . The naive pro¬ 
posal would be to act only on some subregion R x of spins 
near x, but this is clearly unsatisfactory: such an action 
does not commute with the Hamiltonian at the boundary 
of this subregion, and so creates spurious excitations at 
this boundary. Also, it clearly has no projective charac¬ 
ter, since the underlying spins are honest linear represen¬ 
tations of G. Instead, we define U g as follows. 

Our operator U g will be designed to act non-trivially 
only when the state being acted on looks like the ground 
state at the boundary of R x - i.e. its reduced density 
matrix is the same as the ground state’s; if not, we define 
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FIG. 12: Local application of g G G for a trivial 
permutation of the anyon types. The initial system 
contains two anyons of types a and a. A ^-defect and 
g~ x defect are nucleated, braided around the a anyon 
and annihilated, leaving a branch cut surrounding the 
anyon. The symmetry g G G is then applied to all 
degrees of freedom inside defect’s path, thereby 
resulting in a low energy state of the original 
Hamiltonian. 


Ug to act as 0. When the state does look like the ground 
state at the boundary of R x , we act by braiding a twist 
defect of g around this boundary and acting with U 9 
inside R x . To explain this, let us first explain what we 
mean by a twist defect. 

A twist defect is a modification of the Hamiltonian 
along a branch cut terminating at the defect. The lo¬ 
cal terms in the Hamiltonian which straddle the branch 
cut then get twisted by g , as illustrated in Fig [l3j For 
example, if G = Z2 is the Ising symmetry, a term crfcrj 
coupling two sites i and j on opposite sides of the branch 
cut flips sign: crfcrj —>• — crfcrj. Note that this procedure 
is ambiguous near the endpoints of the cut, where the 
“sides” of the cut are not well defined - such microscopic 
choices turn out to lead to discrete gauge ambiguities in 
the definitions below. 


Now to define Ug, we simply nucleate a pair consisting 
of a ^-defect and a g~ x defect, connected by a branch cut, 
and then grow the branch cut to encircle the region R x 
(Fig. 12). This process defines a continuous path in the 
space of gapped Hamiltonians, and Ug acting on a state 
1 ^) is then defined by adiabatically following |^) through 
the branch cut nucleation and growth process, and then 
acting by U g on all of the spins inside the region R x . We 
fix the overall phase of Ug by demanding that Ug |0) = 
|0), where |0) is a state containing no anyons. Because 
the Hamiltonian is gapped, we expect that the adiabatic 
continuation can be approximated exponentially well by 
an operator acting in a thin ribbon (several correlation 
lengths thick) around the boundary of R x . Hence Ug acts 
like U g in the interior of R x , like the identity outside R x , 
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FIG. 13: The Hamiltonian H = JA Hi , where the Hi 
are local terms. The introduction of the ^-defect 
modifies any Hi (blue) which straddles the branch cut 
(green). In such a case, Hi is written as a sum of terms 
which are products of operators acting on the spins to 
the left and right of the branch cut, and then 
conjugates only the right ones by the action of g E G. 
Note that this definition is unambiguous only far away 
from the ^-defect itself. 


and interpolates smoothly between these options within 
a few correlation lengths of the boundary of R x . 

Because we made some arbitrary choices - e.g. the na¬ 
ture of the Hamiltonian at the defect core - we want to 
see to what extent Ug is well defined. A naive argument 
would say that two different such constructions of Ug 
have the same commutation relations with all local oper¬ 
ators near x, and their overall phases are fixed to be the 
same, so they must be equal. This is not quite correct, 
however, because there exist excitations near x - namely 
anyons - which cannot be created with local operators, 
and are sensitive, via braiding, to possible topological 
charge bound to the encircling ^-defect. That is, if we 
had chosen a different nucleation process with different 
defect cores, it is possible that it would differ from the 
first by an additional abelian anyon stuck on the defect, 
which could measure topological charge at x via braiding 
(a non-abelian anyon stuck to the defect would be more 
problematic, and may not lead to a well-defined Ug ; we 
evade this possibility by insisting on a fully gapped de¬ 
fect). Now, there is no way to unambiguously measure 
the abelian anyon charge on a single defect, so we cannot 
specify the abelian anyon charge carried by the defect in 
our procedure - the best we can do is to specify the local 
G-action up to the gauge equivalence: 


Ug\a) ->• e 2ni{a ’ a ^u^\a) (Bl) 

Here | a) is any state with total anyon charge a inside 
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the branch cut loop, a(g) is an arbitrary abelian anyon, 
and the angular brackets denote the full braiding phase of 
a(g) around a (i.e. e 27r d«(^) 5 a > — S a ( g ^ a is their S'-matrix 
element). 


2 . Projective nature of local action of G on anyons 

While the global action of G satisfies the group re¬ 
lation UgUh = Ugh , the corresponding local operators 
Ug Ug h might only satisfy it up to a Berry’s phase, 
which depends on g,h, and the anyon charge a near x. 
To demonstrate this, let us form the operator that would 
measure this Berry’s phase: 


V* h = U x gh {UZ)-\U x g )~ l (B2) 

Now, consider any single site operator Si acting on a spin 
i well inside the local neighborhood R x of x - i.e. inside 
the encircling branch cut. Its commutation relations with 
Ug are the same as its commutation relations with U g - 
this is just the statement that U g acts just like U g inside 
the branch cut region. 



V g h acts only on these spins. 


FIG. 14: Local action of G on anyon a at position x. 
While the operators U gl and Ug h act non-trivially 
in the entire disk, Vg h = L^(L r ^) _1 (C/^) _1 can be 
approximated by an operator that acts only on the 
spins within the narrow surrounding grey ribbon region. 


Si must commute with since the latter acts as 

Ugh{Uh)~ 1 {U g )~ 1 = 1 well inside R x . Thus, since Vg h 
commutes with all observables inside the cut, we can ap¬ 
proximate it by an operator that acts only on spins in 
a thin ribbon region around the cut itself, insofar as its 
action on excitations well inside the cut is concerned - see 
Fig. 14 (the approximation error is exponentially small 
in the thickness of the ribbon divided by the correlation 
length). 

When an interior excitation is created by a local opera¬ 
tor acting on the vacuum, i.e. is of the form | loc) = Y|0), 
with Y a local operator acting only on spins inside the 
cut, we have 


V£ h \loc) = V£ h Y |0> = YV£ h \0 ) = \loc); 

Vg h and Y commute because they are spatially sepa¬ 
rated. Thus there are no Berry’s phases for excitations 
that can be made locally, as expected. However, consider 
now a state | a) which contains an anyon a at x (and no 
other quasiparticles within the branch cut region). | a) 
cannot be nucleated locally from the ground state, but 
rather must be created using a string operator that in¬ 
tersects the branch cut, and hence might fail to commute 
with Vg h , resulting in a Berry’s phase, 


V* h \a) = e i ^’ h ^\a), (B3) 

To compute this action of V gh -> n °f e first that it must 
be by an a-dependent phase, since Y gh must be propor¬ 
tional to the identity within each topological superselec¬ 
tion sector, by the previous argument. Furthermore, if c 
is in the fusion product of a and 6, the phase for c must 
be the sum of the phases for a and 6, because there exists 
a local operator that fuses a and b into c. Mathemati¬ 
cally, any such function from the anyons to U(l) must 
be equivalent to the braiding of an abelian anyon - which 
we call u(g,h) - around the anyon in question. Note that 
Vg h has support on a thin ribbon surrounding the anyon 
in question, so this is certainly an allowed possibility. We 
thus have 


a) = S u{g , hU \a) (B4) 

where S^g^.a is the full braid ^-matrix element. 

Note that there are constraints among the 
u(g,h). Specifically, consider the operator 

W = UfghiUfi)- 1 ^)- 1 ^)- 1 . We can write it 
in the following two ways: 


w — v x v x 

W - V fg,h V f,g 


and 


W = V f,9h U f V 9%( U f)~ 1 - 

The first tells us that W acts by braiding cj(/, g) and then 
braiding cj(fg,h), while the second says that it is equal 
to the composition of braiding operations of u(g, h ) and 
cj(/, gh) (the conjugation by Uj does not affect anything, 
because Uf does not change the topological charge a that 
is being acted on). Thus, we have 

u(fg, h)u(f, g) = u>(f, gh)uj(g, h) 

or 

Those acquainted with group cohomology will recognize 
the above as the cocycle equation 
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doj(f,g,h) = l (B5) 

Furthermore, we can use the definition (|B2 |) o f V^ h to 
see that under the inherent gauge ambiguity ( |Bl| ), u(g, ft) 
must transform as 

u(g, h) w(g, h)a(g)a(h)a(gh)~ 1 


which can also be re-written in group cohomology terms 
as 


u(9, h) ->■ u(g, h ) • da(g, h) 


(B6) 


Functions u(g,h) sati sfyin g eq. ( |B5| ), modulo the gauge 
equivalences in eq. (B6), form a finite abelian group 
called the second group cohomology H 2 (G , *4). Here A is 
the abelian group formed by anyons of quantum dimen¬ 
sion 1. The equivalence class of uj in this set is denoted 
[cj] G H 2 (G,A), and is an invariant which encapsulates 
the projective properties of the anyons in our symmetry 
enriched topological phase. 


Appendix C: Deformations of the defect fusion 
product 



The way that these three equations were obtained from 
the prescription in figur e [10| is illustrated graphically in 
figures [f~5^a) , |l5|fo) , and |16| respectively. 



Here we give a precise description of the deformed de¬ 
fect product defined in figure [lOj Recall that the data 
defining the original defect product of two defects of type 
g and ft is a collection of fusion rules for g and ft supers¬ 
election sectors, and F-matrices involving the triple con¬ 
sisting of a sector of g, a sector of ft, and an anyon, in any 
order. In the following we assume no fusion multiplicity, 
which makes the formulae simpler; the general case is 
not any harder but the formulae are more cumbersome 
to write out. 

Suppose we deform the defect product by a g which 
we denote a for short. Then the deformed fusion rules 
are simply of the form x x y = z where x, y, and z 
are sectors of the g, ft, and gh defects respectively and 
x x y = az is an allowed fusion rule of the undeformed 
defect product. Notice that sometimes, as in the case of 
our Z 4 gauge theory with twisted Z 2 symmetry example 
with a = (0,2), the deformed and undeformed fusion 
rules are identical. To see the difference between the two 
fusion products, we then have to examine the F-matrices. 

Figure [lO] is a shorthand for describing the deformed F 
matrices, which we denote F in contrast to the original 
ones, which we call F. The precise formulae for F in 
terms of F are 


fib,x,y 

X 

zpb,x,y 

X 

rpb,z' ,a 

z 

z'a 

r za 

z' 

r z 


FIG. 15: Definition of deformed F-matrix. 



FIG. 16: Definition of deformed F-matrix. 


Now let us analyze the specific example of the Z 4 gauge 
theory with twisted Z 2 = l,g symmetry. Our claim is 
that deforming the g x g = 1 fusion rule by the anyon 
a = dg^g = (0, 2) (i.e. the strength 2 flux of Z 4 ) forces the 
braiding process of taking a c = (1,0) charge around the 
g defect twice to acquire a minus sign. Note that braiding 
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7-1 X,C,X 

*ba 


X C X 


V 


V 





ba 


x c x 


x c x 


ba 



FIG. 17: Instance of hexagon equation for the original 
defect fusion product. 



FIG. 18: Instance of hexagon equation for the deformed 
defect fusion product. 


this charge around the g defect only once changes it into 
a c = (3,0) charge, and there is no gauge invariant ampli¬ 
tude associated with this process, at the level of extrinsic 
defects (once the theory is gauged, c and c combine into 
a single non-abelian anyon and a single braid maps the 
appropriate fusion space to itself). 

This relation between deforming F-matrices and de¬ 
forming braiding phases of an anyon around a defect fol¬ 
lows from the hexagon equation illustrated in figures [T7| 
and 18 To see this, first note that all of the F and R 
matrices are phases, as the anyons a, 6 , c are all abelian. 


The equations in figures n7\ and 18 read: 


R x,c 

T-:X,C,X 

ba 

[Fir) 


1 pbac,c 

(C4) 




— 1 / ~ —\ 

-1 7 


px,c 

^iX,C,X 

= h b 

(K x ’ x ) 


i R bc ’ c 

(C5) 


Now using equations |C 1 |C 2 [ and |C3| and the fact that 
the F-matrices for the Z 4 gauge theory are all equal to 1 
we obtain 


-frc,X,X _ 771 c,x,x 

b b — b ba 

px,x,c _ px,x,c ^Jp,a 

-frx,c,x rix,c,x 
b b ~ *bci 


so that 


J^X, c J^X, c 

Rx,cRx,c 


R c,a R bc,c (R^ ac ’ c ) 


(C 6 ) 

(C7) 

(C 8 ) 


(C9) 


This is equal to R c ^ a R a ^ c = — 1 . Using the fact that 
in our convention, braiding in front of the defect gives a 
trivial phase of 1 , this shows that the quantity in equa¬ 
tion |80] acquires a —1 when the defect fusion rules are 
deformed by a = ( 0 , 2 ), as desired. 


Appendix D: Hp(G,H) and group extensions 

Symmetry enriched phases are classified by the group 
extensions of H by G. Such a group extension satisfies 

E/H = G (Dl) 

which means that H appears as a normal subgroup of E. 
Each element of E can be labeled by a pair of elements 
h G H and g G G. However, we note that (e, gi) • (e, < 72 ) 7^ 
( e ,<7i<72) in general, as only the cosets E/H multiply as 
G. Thus, a function uj : G 2 H should be introduced 
such that 


(e,3i) • (e,S2) = (w(ffi,S2),ffiS2) (D2) 

G also has a nontrivial effect on i7, permuting the ele¬ 
ments inside the normal subgroup. This is an automor¬ 
phism of i7, and so we can summarize this effect by 

(e,g) ■ (h,e) = (h p{9 \g) (D3) 

p e Hom(G, Aut(iJ)) (D4) 

Combining these data, the total multiplication rule is 

(hi,gi) ■ (h 2 ,g 2 ) = Oi/i2 (9l) w(fli, 92 ), 9152 ) (D5) 

Associativity imposes a constraint on cj, namely 

w(52,53) p(si) w(5i, 5253) = w(pi, 52)^(3152,53) (D6) 

This is precisely the 2-cocycle condition from a group co¬ 
homology of G, Hp(G,E[). This information will be en¬ 
coded into the plaquette terms that live entirely within 
vertical links in our Hamiltonian. Different cocycles cor¬ 
respond to different group extensions that we can con¬ 
struct, producing distinct SET phases. 
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Appendix E: Explicit derivation of plaquette terms 
after gauging 


where we recall that a' L \rj' L ) = M(rj' L )\rj' L ). Focusing on 
ftp(g) we observe that we can simplify it as follows 


As it appears often in this derivation, let gi = nusr. 
Introducing our operator following our prescription in 
eq. [48j we find the plaquette terms are modified as 

Tr# JJ ai{g) -A Tt h ,g JJ 
LeP LeP 

= Tr H II di (flife) - 1 S') ( } TroJJai 

i= 1 i=l 

= Tr ff Op( 5 )n(o') pfe - l) T^lJa* 

2=1 2=1 

=if (&> fe ) -1 £) p(9 


j 

ttp(g) = JR 


pRi-l) 


2 = 1 


(&,(&) 1 < 7 ) 

(ft, (ft) -1 5 ) 

J 

=(9i~i,9i)=tip 


U (gi-i,gi) 


2=1 


which is manifestly independent of g , and thus greatly 
simplifies our plaquette term. 
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